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PREFACE 


The  work  described  in  this  report  was  authorized  under  Project  No.  206023. 
Low-Level  Toxicology.  The  work  was  started  in  June  2003  and  was  concluded  in  January  of 
2004.  The  experimental  data  is  contained  in  laboratory  notebook  02-0109  and  on  compact  discs. 
Raw  data  and  the  final  report  from  this  study  will  be  stored  in  the  Toxicology  Archives.  Building 
E-3150,  Aberdeen  Proving  Ground,  MD. 

In  conducting  this  study,  investigators  adhered  to  the  "Guide  for  the  Care  and  Use 
of  Laboratory  Animals,”  National  Institutes  of  Health  Publication  No.  86-23,  1985.  as 
promulgated  by  the  committee  on  Revision  of  the  Guide  for  Laboratory  Animal  Facilities  and 
Care  of  the  Institute  of  Laboratory  Animal  Resources,  Commission  of  Life  Sciences,  National 
Research  Council,  Washington,  D.C.  These  investigations  were  also  performed  in  accordance 
with  the  requirements  of  AR  70-18,  "Laboratory  Animals,  Procurement,  Transportation,  Use, 
Care  and  Public  Affairs,”  and  the  U.S.  Army  Edgewood  Chemical  Biological  Center  Institutional 
Animal  Care  and  Use  Committee  (IACUC),  which  oversees  the  use  of  laboratory  animals.  This 
project's  assigned  IACUC  protocol  No.  02-341  was  approved  on  6  August  2002. 

Use  of  either  trade  names  or  manufacturers'  names  in  this  report  does  not 
constitute  an  official  endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for 
purposes  of  advertisement. 

This  report  has  been  approved  for  public  release.  Registered  users  should  request 
additional  copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should 
direct  such  requests  to  the  National  Technical  Information  Service. 
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ESTIMATING  LETHAL  AND  SEVERE  TOXIC  EFFECTS  IN  MINIPIGS 
FOLLOWING  10,  60,  AND  180  MINUTES  OF  WHOLE-BODY  GB  VAPOR  EXPOSURE 


1.  INTRODUCTION 

Throughout  history,  research  studies  on  nerve  agents,  have  been  conducted  on  all 
possible  routes  of  human  exposure:  intravenous,  percutaneous,  intramuscular,  subcutaneous, 
intragastric,  and  inhalation.  While  these  studies  provide  valuable  basic  information  on  nerve 
agent  intoxication,  accurate  extrapolation  from  one  route  of  entry  to  another  is  difficult. 
Additionally,  the  progression  of  toxic  signs  via  one  exposure  route  may  be  different  from 
another.  For  instance,  in  an  inhalation  exposure  of  sarin  (GB)  vapor,  constriction  of  the  pupils  is 
the  first  noticeable  sign  of  exposure  while  in  GB  exposures  by  the  subcutaneous  or  intravenous 
routes,  pupil  constriction  occurs  well  after  the  other  signs  of  exposure  or  occurs  irregularly. 

The  likely  routes  of  exposure  for  the  warfighter  in  a  battlefield  situation  would  be  limited  to 
percutaneous  and  inhalation.  Percutaneous  exposure  is  possible  through  contact  with 
contaminated  equipment  or  aerosol  created  from  an  attack.  However,  in  such  an  attack, 
inhalation  of  nerve  agent  vapor  would  be  the  primary  route  of  exposure.  If  GB  vapor  were  used, 
the  chances  of  a  percutaneous  exposure  are  slim.  GB  has  a  vapor  pressure  of  2.9  mm  Hg  at 
25  °C  and  is  likely  to  evaporate  long  before  it  can  be  absorbed  through  the  skin.  Indeed,  open  air 
testing  has  demonstrated  that  GB  has  minimal  effectiveness  in  humans  via  the  percutaneous 
route  of  exposure.2  Therefore,  the  primary  hazard  posed  by  GB  is  through  inhalation. 

Unfortunately,  the  safety  and  logistics  of  performing  an  inhalation  study  have 
been  limiting  factors  on  the  type  of  data  collected.  Historically,  data  collected  during  an 
exposure  have  been  subjective  and  limited  to  clinical  observations.  Additionally,  most  ol  the 
data  belonged  to  the  post-exposure  time  period;  therefore,  physiological  changes  occurring 
during  a  whole-body  inhalation  exposure  have  been  insufficiently  documented.  The  methods 
described  by  Hulet  et  al.  have  made  it  possible  to  collect  data  during  a  real  time  inhalation  GB 
exposure.3 


Traditionally,  the  military  and  other  organizations  dealing  with  inhalation 
toxicology  have  accepted  Haber's  principle  of  dosage,  the  product  of  concentration  (C)  and  time 
(T),  as  constant  over  time  when  assessing  the  impact  of  nerve  agent  vapor  exposures.  Haber  s 
rule  was  used  to  extrapolate  dose/response  data  (based  upon  relatively  short  exposure  times)  to 
predict  response  probabilities  involving  longer  exposure  times.  However,  this  concept  is  now 
considered  inadequate  for  assessing  the  biological  effects  of  exposure  to  many  acutely  toxic 
gases  and  aerosols.5  Recent  efforts3  6  7  have  resulted  in  data  including  low  concentration 
exposures  over  long  periods,  which  can  best  be  described  with  a  toxic-load  model.'  For  even  a 
clear  toxicological  endpoint  as  lethality,  historical  assumptions  used  to  extend  the  prediction  of 
exposures  over  time  have  been  shown  to  be  overly  conservative  for  GB,  which  is  the  best-studied 
agent. 


In  the  toxic-load  model,  the  median  effective  dosage  (ECT50)  increases  with 
exposure  time  in  a  non-linear  relationship  and  the  data  can  be  fit  to  the  toxic  load  equation. 
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CnT  =  k. 


C  is  the  concentration, 

T  is  the  exposure  duration,  and 

n  is  the  toxic  load  exponent,  which  is  dependent  on  the  vapor  or  exposure  scenario. 

To  develop  models  for  predicting  the  probability  of  toxicity  from  low-level  nerve 
agent  exposures  for  different  concentrations  and  durations  of  exposure,  additional  data  from  a 
non-rodent  species  are  needed.  Pigs  have  been  found  to  be  similar  in  anatomy  and  physiology  to 
humans.8  The  intent  of  the  current  studies  is  to  estimate  the  lethal  concentrations  of  the  nerve 
agent  sarin  (GB)  as  a  function  of  exposure-duration  in  the  Gottingen  minipig.  These  results 
helped  determine  whether  LCT50  and  ECT50  (severe)  values  change  with  time,  in  the  minipigs  as 
was  found  in  previous  studies  with  rats.Q 


2.  MATERIALS  AND  METHODS 

2.1  Gottingen  Minipigs. 

Male  and  female  Gottingen  minipigs  were  obtained  from  Marshall  Farms  USA 
(North  Rose,  NY).  Upon  arrival  at  the  testing  facility,  the  minipigs  underwent  an  initial  health 
examination  by  the  attending  veterinary  staff.  The  pigs  were  then  quarantined  for  at  least  three 
days.  After  this  time  the  involved  research  personnel  familiarized  the  pigs  to  various  procedures 
that  included  daily  handling,  change  of  location  within  the  animal  facilities  and  adaptation  to  a 
sling  apparatus.  While  the  animals  were  in  their  cages  their  existence  was  enriched  by  human 
interaction  and  unfettered  access  to  play  toys  (hanging  chains,  bunny  balls)  or  food  treats  on  a 

daily,  rotating  schedule. 

2.2  Surgical  Procedure. 

A  more  thorough  description  of  surgical  procedures  can  be  found  in  Hulet  et  al, 

2006. 3  Briefly,  the  surgical  site  (lateral  neck  from  mandible  to  shoulder  and  mid  dorsally 

between  the  shoulder  blades  of  a  minipig)  was  prepared  for  aseptic  surgery  by  close-clipping  the 

area  and  applying  a  surgical  scrub  (chlorhexidine  or  povidone-iodine)  followed  either  by  an 

application  of  isopropyl  alcohol  or  a  sporicidal  agent.  The  surgical  area  was  then  covered  with 

sterile  drapes  and  the  minipig  positioned  for  surgery  on  a  heated  surgical  table.  The  animal  was 

evaluated  throughout  the  procedure  using  an  EKG  monitor,  temperature  probe,  pulse  oximeter  # 

and  respiratory  monitor. 

m 

A  silicone  catheter  (Bard  access  systems.  6.6  or  9.6  Fr.),  impregnated  with 
heparin  and  antimicrobial  agent,  was  implanted  in  an  external  jugular  vein  and  advanced  to  the 
anterior  vena  cava  or  right  atrium.  A  subcutaneous  tunnel,  extending  from  the  surgical  site 
adjacent  to  the  jugular  vein  to  the  exit  site  in  the  dorsal  midline,  was  created  with  a  hollow 
stainless-steel  rod.  The  catheter  was  filled  with  sterile  heparin  saline  (1%),  grasped  and  pulled 
through  from  the  dorsum  to  the  ventral  neck  incision  with  at  least  6  in.  of  catheter  remaining 
external  to  the  skin.  The  catheter  position  was  adjusted  so  that  blood  samples  could  be  readily 
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obtained.  The  catheter  was  secured  by  tying  at  least  2  sutures  around  the  vein.  A  loop  of 
catheter  leading  from  the  vein  was  also  secured  to  the  subcutaneous  tissues  using  sutures.  Once 
the  catheter  was  appropriately  adjusted,  it  was  secured  at  the  dorsal  exit  site  and  the  incisions 
closed.  The  catheter  was  locked  with  1%  sterile  heparin  saline.  Antibiotic  ointment  was  placed 
on  both  incisions.  Postoperatively,  the  minipig  was  given  analgesics  (buprenorphine  0.01  - 
0.05mg/kg.  BD)  for  at  least  24  hr  and  subsequently,  if  indicated. 

The  minipigs  were  allowed  at  least  3  days  for  recovery  from  the  surgical 
implantation  of  the  indwelling  catheters  before  being  used  for  exposure  to  nerve  agent  vapor. 
During  the  waiting  period,  the  vascular  access  ports  on  the  pigs  were  flushed  with  heparinized 
saline,  as  needed.  During  the  agent  exposures,  the  catheters  were  maintained  by  a  continuous 
intravenous  infusion  of  lactated  Ringers  solution. 

2.3  Blood  Sample  Collection. 

“Real  time”  blood  samples  were  drawn  via  the  indwelling  jugular  catheters  to 
assess  cholinesterase  inhibition  and  GB  uptake  via  regeneration  assays.10  Blood  samples  were 
taken  just  prior  to  exposure  and  periodically  throughout:  approximately  every  2  min  during  the 
10-min  exposure,  every  15  min  during  the  60-min  exposure,  and  every  20  min  during  the  180- 
min  exposure.  The  total  volume  of  blood  drawn  did  not  exceed  1%  of  the  body  weight  of  a 
minipig  over  a  1-week  span.  The  drawn  sample  volumes  were  replaced  by  an  equivalent  volume 
of  Lactated  Ringers. 

Assays  for  AChE  and  BChE  activity  were  performed  on  whole-blood  samples. 
Ten  uL  of  whole  blood  was  added  to  a  disposable  borosilicate  glass  tube  (Chase  Scientific  Glass, 
Rockwood,  TN)  containing  2000  uL  of  distilled  water.  Two  hundred  uL  of  0.69  mM  phosphate 
buffer  at  pH  7.4  (EQM  Research,  Cincinnati,  OH)  was  then  added  to  each  tube.  The  tubes  were 
vortexed  and  allowed  to  sit  at  room  temperature  for  20  min.  Two  hundred  uL  of  the  sample 
solution  from  each  tube  was  transferred  to  individual  wells  on  a  96-well  plate.  Twenty-five  uL 
of  30  mM  5,  5-dithiobis-2-nitrobenzoic  acid  (DTNB)  was  added  to  each  well.  The  plate  was 
covered,  and  incubated  at  37  °C  for  15  min. 

For  the  determination  of  AChE  activity,  25  uL  of  a  solution  containing  10-mM 
acetylthiocholine  and  200-uM10-(a-diethylaminopropionyl)-phenothiazine,  a  specific  inhibitor 
of  butyrylcholinesterase  (EQM  Research,  Cincinnati,  OH),  was  added  to  the  appropriate  wells  of 
the  96-well  plate.  For  the  determination  of  BChE  activity,  25  uL  of  a  solution  containing  20-mM 
butyrylthiocholine  (EQM  Research,  Cincinnati.  OH)  was  added  to  the  appropriate  wells  of  the 
96-well  plate.  The  plate  was  shaken  briefly  to  ensure  mixing  of  the  reagents  and  read  at  450  nm 
and  37  °C  using  a  SpectraMax  Plus  microplate  spectrophotometer  (Molecular  Devices  Corp., 
Sunnyvale,  CA)  for  10  min  and  analyzed  using  SoftMax  Pro  LS  version  4.3  software.  AChE  and 
BChE  activity  values  were  expressed  as  units  of  activity  per  mL  of  whole  blood  (U/mL). 

2.4  Inhalation  Chamber. 

The  minipig  whole-body  exposures  were  conducted  in  a  1000-L  dynamic  airflow 
inhalation  chamber.  The  Rochester  style  chamber  was  constructed  of  stainless  steel  with  glass  or 
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Plexiglas  windows  on  each  of  its  6  sides.  The  interior  of  the  exposure  chamber  was  maintained 
under  negative  pressure  (0.25-0.30"  FTO),  which  was  monitored  with  a  calibrated  magnehelix 
(Dwyer,  Michigan  City,  IN).  A  thermoanemometer  (Model  8565,  Alnor,  Skokie,  IL)  was  used  to 
monitor  chamber  airflow  at  the  outlet. 

Two  sampling  methods  were  used  to  monitor  and  analyze  the  GB  vapor 
concentration  in  the  exposure  chamber.  The  first  method  was  a  quantitative  technique  using 
solid  sorbent  tubes  (Tenax/TA)  to  trap  GB,  followed  by  thermal  desorption.  The  second  method 
was  a  continuous  monitoring  technique  using  a  phosphorus  monitor  (HYFED.  Model  PA260  or 
PH262,  Columbia  Scientific,  Austin,  Texas).  Output  from  the  HYFED  monitor  provided  a 
continuous  strip  chart  record  of  the  rise,  equilibrium,  and  decay  of  the  chamber  vapor 
concentration  during  an  exposure. 

All  air  samples  were  drawn  from  the  middle  of  the  chamber  and  solid  sorbent  tube 
samples  were  drawn  after  the  chamber  attained  equilibrium  (t99).  The  HYFED  monitored  the  GB 
vapor  concentration  during  the  entire  run.  Solid  sorbent  tube  samples  were  drawn  from  the 
chamber  approximately  every  10  min  with  each  sample  draw  lasting  1  to  5  min  depending  upon 
chamber  concentration  and  duration  of  exposure.  All  sample  flow  rates  for  the  solid  sorbent  tube 
systems  were  controlled  with  calibrated  mass  flow  controllers  (Matheson  Gas  Products, 
Montgomeryville,  PA).  Flow  rates  were  verified  before  and  after  sampling  by  temporarily 
connecting  a  calibrated  flow  meter  (DryCal®,  Bios  International,  Pompton  Plains,  NJ)  in-line  to 
the  sample  stream.  Physical  parameters  (chamber  airflow,  chamber  temperature,  and  relative 
humidity)  were  monitored  during  exposure  and  recorded  periodically. 

2.5  Solid  Sorbent  Tube  System. 

The  automated  solid  sorbent  tube  sampling  system  consisted  of  four  parts: 

(1)  A  heated  sample  transfer  line, 

(2)  A  heated  external  switching  valve, 

(3)  A  thermal  desorption  unit,  and 

(4)  A  gas  chromatograph  (GC). 

A  steel  sample  line  (1/16"  o.d.  x  0.004"  i.d.  x  6'  length)  extended  from  the  middle 
of  the  chamber  to  an  external  sample  valve.  The  sample  line  was  commercially  treated  with  a 
silica  coating  (Silicasteel®  Restek.  Bellefonte,  PA)  and  covered  with  the  heated  (60  °C)  sample 
transfer  line  (CMS,  Birmingham,  AL).  The  combination  line  coating  and  heating  minimized 
nerve  agent  absorption  onto  the  interior  surfaces  of  the  sample  line.  From  the  transfer  line,  the 
sample  entered  a  heated  (125  °C)  6-port  gas-switching  valve  (UWP,  Valeo  Instruments, 

Houston,  TX).  In  the  by-pass  mode,  vapor  from  the  chamber  continuously  purged  through  the 
sample  line  and  out  to  a  charcoal  filter.  In  the  sample  mode,  the  gas  sample  valve  redirected 
nerve  agent  vapors  from  the  sample  line  to  a  Tenax  TA/Haysep  sorbent  tube  (60-80  mesh) 
located  in  the  thermal  desorption  unit  (ACEM-900,  Dynatherm  Analytical  Instruments,  Kelton. 
PA).  Temperature  and  flow  programming  within  the  Dynatherm  desorbed  nerve  agents  from  the 
sorbent  tube  directly  onto  the  GC  column  (RTX-5,  30  m,  0.32  mm  i.d.,  1-mm  thickness),  which 
was  then  followed  by  flame  photometric  detection  (FPD  -  phosphorus  mode). 
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The  solid  sorbent  tube  sampling  system  was  calibrated  by  the  direct  injection  of 
external  standards  (GB  pg/ml)  into  the  heated  sample  line  of  the  Dynatherm.  This  way,  injected 
nerve  agent  standards  were  put  through  the  same  sampling  and  analysis  stream  as  the  chamber 
samples.  A  linear  regression  fit  (r2  =  0.999)  of  the  standard  data  was  used  to  compute  the  GB 
concentration  of  each  chamber  sample. 

2.6  Chemicals. 

Isopropyl  methyl  phosphonofluoridate  (GB)  was  used  for  all  the  vapor  exposures 
in  this  study.  Chemical  agent  standard  analytical  reagent  material  (CASARM)-grade  GB  2035 
•  (lot  #  GB-U-6814-CTF-N)  was  verified  (usually  98.3  +  0.48  wt.  %  pure  as  determined  by 

quantitative  3IP-NMR)  and  stored  in  sealed  ampoules  containing  nitrogen.  The  ampoules  were 
opened  as  needed  either  to  prepare  external  standards  or  to  be  used  as  neat  agent  for  vapor 
generation.  All  external  standards  for  GB  vapor  quantification  were  prepared  on  a  daily  basis. 
Triethylphosphate  (99.9%  purity),  obtained  from  Aldrich  Chemicals,  Milwaukee,  WI,  was  used 
as  the  internal  standard  for  the  GB  purity  assays. 

Analysis  for  agent  impurities  was  conducted  using  acid-base  titration  as  well  as 
Gas  Chromatography/Mass  spectrometry  (GC-MS)  and  NMR.  Acid-base  titration  has  been 
found  to  show  the  following  impurity  percentages  based  on  mole  ratios: 

GB  ANALYSIS 


Compound  Mole  %  Calculated  Wt  % 

Methylphosphonofluoridic  acid  (Fluor  Acid)  0.3  0.2 

Diisopropyl  methylphosphonate  (DIMP)  0.2  0.3 

Methylphosphonic  difluoride  (DF)  0.2  0.2 

GC-MS  positively  identified  DIMP,  Diisopropyl  phosphonofluoridate,  Tributylamine,  and 
Isopropyl  ethylphosphonofluoridate,  but  did  not  quantify  the  amounts.  Tributylamine  was  also 
confirmed  using  *H  NMR  with  a  concentration  of  <  0. 1  weight%  of  GB. 

2.7  Vapor  Generation. 

Saturated  GB  vapor  streams  were  generated  by  flowing  nitrogen  carrier  gas 
through  a  glass  vessel  (multi-pass  saturator  cell)  that  contained  liquid  GB.  The  saturator  cell 
consists  of  a  100  mm  long,  25-mm  o.d.  cylindrical  glass  tube  with  two  (inlet,  outlet)  vertical 
7-mm  o.d.  tubes  connected  at  each  end.  The  main  body  of  the  saturator  cell  contains  a  hollow 
ceramic  cylinder  that  serves  to  increase  the  contact  area  between  the  liquid  nerve  agent  and  the 
nitrogen.  The  saturator  cell  allows  nitrogen  to  make  three  passes  along  the  surface  of  the  wetted 
ceramic  cylinder  before  exiting  the  outlet  arm  of  the  glass  cell.  The  saturator  cell  body  was 
immersed  in  a  constant  temperature  bath  so  that  a  combination  of  nitrogen  flow  and  temperature 
could  regulate  the  amount  of  nerve  agent  vapor  going  into  the  inhalation  chamber.  The  entire 
apparatus  was  contained  within  a  generator  box  that  was  mounted  at  the  top  of  the  inhalation 
chamber.  Typically,  the  saturator  cell  was  loaded  with  2  to  4  ml  of  liquid  nerve  agent 
(CASARM  grade).  To  maintain  the  integrity  of  the  liquid  nerve  agent  within  the  cell,  a 
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continuous  low  flow  rate  (1  to  2  ml/min)  nitrogen  stream  was  used.  This  setup  was  capable  of 
precisely  generating  GB  vapor  over  a  concentration  range  of  0.001-2.0  mg/m3. 

2.8  Sling  Apparatus. 

A  sling  was  used  to  restrain  each  minipig  during  the  exposure  to  the  GB  nerve 
agent  vapor.  The  frame  of  the  sling  was  constructed  of  airtight  stainless  steel  pipe  and 
Swagelok™  fittings.  The  slings  were  custom  designed  (Lomir  Biomedical,  Inc.,  Malone.  NY  or 
Canvas  and  Awning  supplies.  White  Marsh,  MD)  to  fit  the  build  and  size  of  the  minipigs.  The 
body  of  each  sling  was  made  of  canvas,  which  contained  2  holes  on  each  side  for  legs  so  that  it 
fitted  comfortably  around  the  pig.  The  sling  also  had  two  2  straps  that  secured  over  the 
shoulders  and  hips.  A  muzzle  harness  was  placed  over  the  snout  and  secured  both  laterally  and 
ventrally  to  the  stainless-steel  framing  in  order  to  prevent  the  pig  from  moving  its  head  from 
side-to-side.  The  harness  was  fitted  so  that  it  did  not  interfere  with  the  pig's  ability  to  breathe. 

2.9  Design  and  Data  Analysis. 

To  determine  the  progression  of  experimental  exposure  concentrations,  the  up- 
and-down  method  with  an  assumed  probit  slope  of  10  was  used1 '.  For  this  study,  the  binary 
response  used  for  executing  this  method  was  dependent  on  the  survival  of  the  minipig  for  24  hr 
after  exposure  to  the  nerve  agent.  The  signs  of  nerve  agent  exposure  were  designated  as 
moderate,  severe,  or  lethal.  A  minipig  was  classified  as  having  severe  signs  of  exposure  if  it 
were  gasping,  prostrated,  collapsed,  or  convulsing.  Muscle  tremors,  salivation,  lacrimation,  or 
miosis  constituted  a  moderate  exposure. 

The  method  of  maximum  likelihood  estimation  (MLE)  was  used  on  the  resulting 
quantal-response  data  to  calculate  LCT50  values  (and  associated  asymptotic  95%  confidence 
intervals)  and  ECT50  values  for  severe  signs  for  each  of  the  6-gender  exposure-  duration 
groups12.  An  example  of  an  MLE  calculation  is  presented  in  Appendix  A  of  ECBC-TR-450.3 

Up-and-down  experiments  normally  use  6  to  10  subjects,  definitely  not  enough  to 
permit  a  reliable  estimation  of  the  probit  slope.  However,  data  from  several  up-and-down 
experiments  can  be  combined  to  obtain  enough  animals  (at  least  30)  to  estimate  the  probit  slope. 

The  resulting  dataset  can  then  be  analyzed  via  traditional  probit  analysis  or  ordinal  logistic 
regression  to  obtain  a  probit  slope  estimate.7 13 14  15 

Equations  1  and  2  were  used  to  model  the  response  distribution: 

Yn  =  (Yp  ~  5)  =  *0  +  Mlogio  0+M*°gio  T)+k,  (Gender )+kTS  (logl0  T)(Gender )  (1) 

Yn  =  (Yr~  5)  =  k0  +  kc(\ogl0C)  +  kr(Time)  +  ks  (Gender)  +  kTS  (Time)(Gender)  (2) 

where  Yn  is  a  normit;  Yp  is  a  probit;  the  k's  are  fitted  coefficients;  C  is  vapor  concentration;  and 
both  T  and  Time  represent  exposure-duration.  In  eq  1 ,  exposure-duration  is  treated  as  a 
covariate  (T),  whereas  in  eq  2,  exposure-duration  is  treated  as  a  3-level  factor  (Time).  The 
constant,  krs,  has  6  values,  one  for  each  Time-Gender  combination.  The  constants,  kc  and  kj. 
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are  the  probit  slopes  for  concentration  and  Time,  respectively.  The  toxic  load  exponent,  n,  is  the 
ratio  kc/kj.  If  this  ratio  is  not  different  (with  statistical  significance)  from  1,  then  Haber's  rule  is 
appropriate  for  modeling  the  toxicity.  Otherwise,  the  toxic  load  model  (CnT)  is  the  proper 
approach,  assuming  that  there  is  no  significant  curvature  in  the  data  used  to  fit  the  model. 

Should  significant  curvature  exist,  the  toxic  load  model  is  not  appropriate,  but  it  is  still  superior 
to  Haber's  rule  in  modeling  the  data. 

The  present  protocol  has  exposure  durations  of  10,  60,  and  1 80  min.  For  each  of 
the  exposure  durations,  6  or  7  minipigs  of  each  gender  were  used.  Statistical  analysis  routines, 
contained  within  Minitab®  versions  13  and  14  (Minitab,  Inc.,  State  College  PA),  and  an  in- 
house  developed  spreadsheet  program  were  used  for  the  analysis  of  the  data. 


3.  RESULTS 

3.1  Animals. 

Thirty-eight  pigs  (19  male  and  19  female)  were  exposed  to  concentrations  of  GB 
vapor  to  estimate  LCT50  and  ECT5o  (severe)  values.  An  additional  male  pig  was  used  as  an  air 
control.  At  the  time  of  the  surgeries,  the  20  males  (19  experimental  plus  1  control),  weighed  an 
average  of  10.68  kg  ±  0.26  (SEM)  kg  and  the  19  females,  weighed  an  average  of  10.62  ±  0.21 
(SEM)  kg. 

3.2  Median  Effective  and  Median  Lethal  Dosages. 

The  results  of  the  exposures  were  classified  as  moderate,  severe,  or  lethal  (see 
section  2.9  for  a  description  of  criteria).  The  ordinal  scores  are  listed  in  Table  1.  The  method  of 
maximum  likelihood  estimation  (MLE)  was  used  to  calculate  ECT50  (severe)  and  LCT50  values 
(and  associated  asymptotic  95%  confidence  intervals)  for  each  of  the  6-gender  exposure-duration 
groups.12  MLE  calculations  are  shown  in  Appendix  A.  The  EC50  and  ECT50  (severe)  values,  with 
their  respective  95%  confidence  intervals,  can  be  found  in  Table  2  and  are  plotted  in  Figures  2 
and  4.  The  ECT50  values,  like  LCT50  values,  are  not  constant  over  time.  The  ratio  of  ECT50 
(severe)  values  to  LCT50  values  are  shown  in  Table  3.  The  ratio  of  lethal  to  severe 
concentrations  was  higher  in  female  pigs  (99%  confidence). 

3.3  Statistical  Models  for  the  Probability  of  Lethality. 

To  model  the  probability  of  lethality  as  a  function  of  exposure  concentration, 
exposure-duration,  and  gender,  several  models  were  fit  to  the  quantal  data  shown  in  Table  1. 

The  number  of  pigs  used  for  each  gender  exposure-duration  group  was  not  large  enough  to 
estimate  the  response  distribution.  Instead,  the  response  distribution  was  estimated  using  either 
eq  1  or  2  with  the  data  for  all  38  pigs  (see  Section  2.9).  Ordinal  regression  was  used  to  fit 
various  response  models  (eqs  1  and  2)  to  the  data.  Appendix  B  contains  printouts  of  the 
MINITAB®  results. 
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Pig  63  may  be  considered  an  outlier.  Tables  4  (analyses  without  pig  63)  and  5 
(analyses  with  pig  63)  give  the  probit  slopes  and  toxic  load  exponents  obtained  from  various 
ordinal  regression  model  fits.  The  recommended  best  model  fit  is  model  L5  (without  pig  63): 

Yn  =  constant  +  12.4  log  10(C)  +  9.0  loglO(T)  -  .605  Sex 

where  the  constant  depends  on  the  effect  (severe  or  lethal)  and  Sex  is  coded  as  1  for  male  and  -1 
for  female. 


The  value  of  the  toxic  load  exponent  (n  =  kC  /  kT)  was  essentially  independent  of 
the  model  used.  The  toxic  load  exponent  of  model  L5  was  1.38  with  a  95%  confidence  interval 
of  1.24  to  1.52  when  pig  63  was  excluded  (if  pig  63  was  included,  n  =  1.37  with  95%  confidence 
interval  of  1.20  to  1.53).  Because  this  interval  did  not  overlap  one.  Haber's  rule  was  not 
considered  an  appropriate  time  dependence  model  for  this  dataset.  Potential  curvature  in  the  data 
was  evaluated  by  inserting  a  (logT)"  term  into  the  model.  This  term  was  found  to  be  statistically 
insignificant,  regardless  of  the  inclusion  or  exclusion  of  pig  63  in  the  analysis. 

For  executing  the  up  and  down  method  in  this  study,  the  probit  slope  on 
concentration.  kC,  was  assumed  to  be  10  (see  Section  2.9).  The  probit  slope  of  the  best  model  fit 
(L5)  was  12.4  with  a  95%  range  of  6.2  to  18.6  when  pig  63  was  excluded  (if  pig  63  were 
included,  the  probit  slope  would  be  9.2  with  a  95%  range  of  4.4  to  14.0).  However,  regardless  of 
the  inclusion  of  pig  63,  all  the  95%  confidence  intervals  for  the  probit  slope  from  the  6  model  fits 
overlapped  10. 

3.4  Gender  Differences. 

The  models  were  tested  for  possible  gender  effects  and  Sex  was  found  to  be  a 
significant  term  (p  =  0.014  in  model  L2  and  p  =  0.013  in  model  L5).  When  pig  63  was  excluded 
from  the  analysis,  male  minipigs  were  noticeably  more  sensitive  than  females.  Gender  was  not  a 
statistically  significant  term  (p  =  0.063  in  Model  L2  and  p  =  0.067  in  Model  L5)  when  pig  63 
was  included  in  the  analysis.  The  gender  term  was  not  statistically  significant  for  any  of  the 
models  where  exposure-duration  was  treated  as  a  covariate.  The  interaction  of  Sex  with  Time 
(Model  LI)  or  Sex  with  logT  (Model  L4)  was  not  statistically  significant,  regardless  of  the 
inclusion  of  pig  63.  The  failure  to  find  statistically  significant  differences  between  the 
interactions  of  gender  and  exposure-duration  may  have  been  due  to  the  low  sample  size 
(n  =  6  to  7).  In  contrast,  when  gender  was  considered,  regardless  of  exposure-duration,  the 
sample  size  was  much  larger  (n  =  38). 

3.5  Baseline  AChE  and  BChE  Activity. 

Baseline  activities  of  red  blood  cell  acetylcholinesterase  (AChE)  and  plasma 
butylcholinesterase  (BChE)  were  assessed  prior  to  nerve  agent  exposures.  Of  the  38  pigs  that 
were  exposed  to  GB  (19  males  and  19  females),  blood  samples  could  not  be  collected  during 
exposure  from  3  males  and  2  females  because  their  catheters  lost  patency  between  the  day  of 
surgery  and  the  day  of  the  experiment.  However,  the  baseline  activities  of  the  male  control  pig 
were  included  for  a  total  of  17  males  and  17  females.  Baseline  AChE  activities  showed  no 
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significant  differences  (p  =  0.286)  when  comparing  male  and  female  minipigs.  However,  female 
minipigs  showed  significantly  less  (p  =  0.022)  baseline  BChE  activity  as  compared  to  male 
minipigs.  In  order  to  increase  the  total  number  of  minipigs  used  for  the  baseline  measurements, 
the  data  from  the  pigs  used  in  the  current  studies  were  combined  with  baseline  measurements 
taken  from  pigs  in  other  ongoing  studies  to  provide  a  total  of  37  females  and  44  males.3  The 
significant  difference  between  baseline  male  and  female  BChE  activity  was  more  pronounced 
(p  =  0.004)  while  there  was  still  no  difference  in  baseline  AChE  activity  (p  =  0.68 1 ). 

3.6  Depression  of  AChE  and  BChE  Activity. 

Depression  of  cholinesterase  (AChE  &  BChE)  activity  was  assessed  during  the 
GB  exposures  by  collecting  blood  specimens  through  the  jugular  catheter.  Of  the  38  pigs  that 
were  exposed  to  GB  (19  males  and  19  females)  blood  samples  could  not  be  collected  from 
3  males  and  2  females  due  to  the  loss  of  patency  in  the  catheter  between  the  day  of  surgery  and 
the  day  of  the  experiment.  For  male  pig  no.  58,  samples  could  only  be  collected  after  exposure 
had  been  concluded  and  the  pig  had  been  removed  from  the  exposure  chamber.  Table  6 
identifies  the  lowest  AChE  and  BChE  values  during  exposure  and  gives  the  absolute  lowest 
values  (during  or  after  exposure)  for  each  pig.  AChE  values  were  decreased  to  8%  or  less  of  the 
baseline  values  before  the  conclusion  of  the  exposures  in  31  out  of  32  pigs.  There  was  very  little 
subsequent  depression  in  AChE  activity  after  the  conclusion  of  the  exposures.  In  fact,  AChE 
depression  in  22  of  32  pigs  had  reached  the  absolute  lowest  value  during  the  course  of  the 
exposure.  Depression  of  BChE  was  variable  but  dropped  below  50%  of  baseline  in  only  one  of 
the  pigs. 

3.7  Rate  of  AChE  Depression  and  GB  Uptake. 

Only  the  blood  samples  collected  during  the  exposure  time  were  utilized  to 
calculate  the  rate  of  AChE  depression.  The  theoretical  dosage  that  each  pig  was  exposed  to,  up 
to  the  point  that  each  blood  specimen  was  drawn,  could  be  calculated.  The  depression  of  AChE 
(expressed  as  a  percent  of  baseline  measurements)  was  then  plotted  versus  exposure  dosage 
(mg.min/m3).  The  depression  of  AChE  activity  versus  dosage  was  best  modeled  using  a 
polynomial  fit.  The  associated  x,  x2  and  r2  values  from  the  polynomial  curve  fits  are  found  in 
Table  7.  The  rate  of  AChE  depression  for  each  pig  at  any  exposure  dosage  could  be  ascertained 
by  taking  the  instantaneous  slope  of  the  curve  at  any  dosage  along  the  X-axis.  There  were  no 
significant  statistical  differences  found  between  the  rates  of  AChE  depression  in  pigs  that  lived 
versus  the  pigs  that  died.  This  was  the  case  whether  all  the  pigs  were  considered  regardless  of 
exposure-duration  or  whether  the  groups  were  broken  down  into  separate  exposure-durations. 

Using  the  same  methods  as  described  above,  multiple  pair  wise  comparison  t-tests 
were  utilized  to  test  for  significant  gender  differences  between  groups  along  a  range  of  dosages 
(Table  8).  When  the  minipigs  were  considered  as  a  group,  regardless  of  the  duration  of 
exposure,  the  rate  of  ChE  depression  was  significantly  different  (p  =  0.043)  between  males  and 
females  at  a  dosage  of  7.0  mg.min/m3  and  slightly  outside  of  statistical  significance  (p  =  0.054) 
at  8.0  mg.min/m3.  At  dosages  lower  than  7.0  mg.min/m3,  the  statistical  differences  became  more 
pronounced  (p  =  0.027  at  Ct  =  5,  p  =  0.012  at  Ct  =  1).  In  all  the  circumstances,  the  statistical 
differences  between  the  rates  of  AChE  depression  in  male  pigs  had  steeper  slopes  than  the 
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female  pigs.  When  the  groups  of  male  and  female  pigs  were  further  broken  down  by  exposure- 
duration,  there  were  no  statistical  differences  at  any  dosage  for  the  180-min  exposures.  There 
were  statistical  differences  between  the  male  and  female  rates  of  AChE  depression  in  the 
60-min  exposures  at  total  dosages  up  to  30  mg.min/m3  (p  =  0.016)  and  for  the  10-min  exposures 
up  to  20  mg.min/m3  (p  =  0.039). 

Only  those  blood  samples  collected  during  the  exposure  time  were  utilized  to 
calculate  the  rate  of  GB  uptake  (see  Section  2.3),  which  was  plotted  versus  exposure  dosage 
(mg/min/m3)  for  each  pig.  These  curves  were  best  modeled  using  linear  fits.  The  associated  x 
(slope)  and  r2  values  for  the  linear  fits  are  found  in  Table  7.  Multiple  pair  wise  comparison 
t-tests  were  utilized  to  test  for  significant  differences  between  groups  (Table  9).  There  were  no 
significant  differences  between  genders  when  males  and  females  were  considered  as  groups 
regardless  of  exposure-duration.  This  was  also  the  case  when  the  groups  were  individually 
broken  down  to  10-,  60-,  or  1 80-min  exposure-durations.  In  contrast,  there  was  a  highly 
significant  difference  (p  =  0.004)  between  the  uptake  rates  of  GB  between  animals  that  survived 
for  24  hr  and  animals  that  eventually  died.  Not  surprisingly,  the  rate  of  uptake  was  higher  in  the 
animals  that  eventually  died. 


4.  DISCUSSION 

4. 1  LCUn  Values. 

The  calculated  LCT50  values  for  pigs  in  the  current  study  are  consistent  with  the 
notion  that  larger  mammals  (pigs,  dogs,  cats,  monkeys)  have  lower  threshold  values  than  smaller 
animals  (mice,  rats,  rabbits).  Larger  mammals  are  also  more  reflective  of  estimated  LCT50 
values  in  humans.  For  instance,  the  10-min  LCT50  values  in  mice,  male  rats  and  rabbits  are  380. 
231  and  1 15  mg.min/m3.6  16  In  comparison,  the  calculated  LCT50S  for  10-min  GB  exposures  in 
monkeys  and  male  cats  are  74  mg.min/m3  and  79  mg.min/m3,  respectively.17 18  Crook  et  al 
calculated  the  LCT50  for  a  10-min  GB  vapor  exposure  in  pigs  to  be  34  mg.min/m3.19  The  LCT50 
for  10-min  exposures  in  male  and  females  pigs  in  the  current  study  are  72.5  and  86.9  mg.min/m3, 
respectively. 


The  approximate  2-fold  difference  in  the  LCT50  values  in  the  two  studies  may  be 
attributed  to  several  variables.  The  Crook  study  used  Yorkshire  pigs  instead  of  Gottingen 
minipigs.  The  age  of  the  Yorkshire  pigs  was  identified  only  by  saying  that  they  were  "just 
recently  weaned."  as  opposed  to  the  sexually  mature  pigs  in  this  study.  The  pigs  in  the  Crook 
studies  were  exposed  in  groups  of  4  or  5  and  were  allowed  free  movement  within  the  chamber, 
while  in  the  current  study,  the  pigs  were  individually  exposed,  restrained,  and  immobile. 
Additionally,  the  methods  for  generating  accurate  concentrations  and  analytically  verifying  those 
concentrations  are  likely  more  reliable  now  than  they  were  50  years  ago.  Despite  the  differences 
in  the  LCT50  values  between  the  current  study  and  the  study  by  Crook,  the  numbers  from  the 
studies  done  on  pigs  are  considerably  closer  to  human  LCT50  estimates  than  the  data  obtained 
from  the  rodent  studies. 
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Bide  et  al.  recently  suggested  an  estimate  for  a  10-min  GB  vapor  exposure  of  57 
mg.min/m3  in  humans.20  This  estimate  was  extrapolated  based  on  data  taken  from  38  historical 
animal  studies  involving  7  species  (none  being  swine),  regardless  of  gender.  The  overview  took 
into  account  the  minute  volume  (MV)  to  body  weight  (B  W)  ratio  for  each  of  the  species.  The 
calculated  MV/BW  ratio  for  humans  was  0.223.  The  next  closest  MV/BW  ratio  of  species  used 
in  the  study  was  0.328  for  dogs.  The  MV/BW  ratio  of  pigs  was  0.225. _l  The  calculated  toxic 
load  exponent  for  the  large  data  set  was  1.38.  Interestingly,  the  calculated  toxic  load  exponent 
for  the  current  study  was  also  1 .38. 

The  LCT50  reported  for  a  10-min  GB  vapor  exposure  in  monkeys  was 
•  74  (62-87  F.l.)  mg.min/m3.17  In  the  Cresthull  study,  monkeys  demonstrating  the  toxicological 

signs  of  collapse  and/or  convulsions  were  categorized  as  “incapacitated"  and  then  an  ICT50 
(incapacitation)  value  for  the  10-min  GB  exposure  was  calculated.  The  ratio  of  the  ICT50  to 
LCT50  was  0.89.  In  the  current  study  the  toxicological  signs  of  collapse  and/or  convulsions 
(along  with  prostration  and  gasping)  were  utilized  to  characterize  a  pig  as  portraying  “severe" 
signs  of  exposure.  The  ratio  of  ECT50  (severe)  to  LCT50  values  in  male  pigs  for  10-,  60-.  and 
1 80-min  exposures  were  0.71,  0.79,  and  0.74,  respectively.  The  ratio  of  ECT50  (severe)  to  LCT50 
for  female  pigs  for  10-,  60-,  and  1 80-min  exposures  were  0.89,  0.89  and  0.84,  respectively. 
Statistically,  the  ratios  of  severe  to  lethal  concentrations  were  higher  in  female  pigs  (99% 
ANOVA  confidence)  than  in  male  pigs. 

The  data  suggest  that  there  is  less  difference  between  severely  toxic  and  lethal 
dosages  in  female  as  compared  to  male  pigs.  While  the  Cresthull  study  used  both  male  and 
female  monkeys,  the  breakdown  of  sexes  only  stated  that  “most  of  them  were  females"  limiting 
the  ability  to  compare  the  two  studies  in  order  to  distinguish  whether  statistical  gender 
differences  exist  in  the  severe  to  lethal  ratios  of  other  species.  The  current  study  is  most  likely 
the  first  to  identify  gender  differences  between  the  ratios  of  severe  to  lethal  effects.  Sommerville 
analyzed  the  data  of  Mioduszewski  et  al.  to  calculate  an  ECT50  (severe)  to  LCT50  ratio  of  0.79 
for  a  vapor  GB  exposure  in  rats.9'15  This  value  was  calculated  for  all  of  the  animals  taken 
together  as  a  single  group,  regardless  of  the  duration  of  exposure  or  gender.  The  rat  data  of 
Mioduszewski  et  al.  is  being  reanalyzed  to  determine  whether  the  data  from  this  study  are 
applicable  to  rats.9 

Crook  et  al.  determined  that  87%  of  the  pigs  that  died  from  the  10-min  GB  vapor 
exposures  did  so  either  during  exposure  or  within  the  first  10  min  after  the  conclusion  of  the 
exposure.19  In  the  current  study,  83%  of  the  pigs  that  died  from  the  10-min  exposure  to  GB  did 
«  so  within  the  first  1 0  min  after  exposure.  Similarly,  72%  of  the  pigs  that  died,  regardless  of 

exposure-duration,  did  so  either  during  the  exposure  or  within  the  first  1 0  min  after  the  exposure. 
These  data  suggest  that  the  toxic  actions  of  GB  occur  because  of  the  inhalation  of  GB  vapor 
rather  than  the  delayed  absorption  of  GB  through  the  skin.  Indeed  it  has  been  demonstrated  in 
open  air  testing  that  GB  has  minimal  effectiveness  in  humans  as  a  nerve  agent  via  the 
percutaneous  route  of  exposure.2 
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4.2 


Gender  Differences. 


In  1998.  a  recommended  change  to  airborne  occupational  exposure  limits 
suggested  that  the  Immediately  Dangerous  to  Life  or  Health  (IDLH)  exposure  guidelines  for 
inhaled  GB  be  lowered."  This  suggestion  was  made  to  correct  the  failure  of  the  existing 
guidelines  to  take  into  account  that  there  may  be  differences  in  sensitivity  to  nerve  agents  based 
on  gender.  The  existing  value  at  the  time  for  a  30-min  exposure  to  GB  of  0.2  mg/m3  was 
lowered  to  0.1  mg/m3.  This  suggestion  was  made  based  on  work  done  by  Callaway  and 
Blackburn  in  which  female  rats  were  found  to  be  as  much  as  twice  as  sensitive  to  the  lethal 
effects  of  inhaled  GB  than  male  rats."  The  significantly  greater  sensitivity  to  inhaled  GB  in 
female  rats  has  subsequently  been  shown  to  occur  over  longer  (240  min)  durations  of  exposure 
as  well.6  The  female  hamster  has  also  been  identified  as  being  more  susceptible  to  GB  vapor 
exposure  than  its  male  counterpart.24  In  contrast,  male  mice  have  been  identified  as  being 
significantly  more  sensitive  than  female  mice  to  GB  vapor  via  inhalation  and  intravenous 
administration."  Given  that  there  are  no  relevant  human  data  available  and  there  is  a  surprising 
lack  of  literature  investigating  gender  differences  in  sensitivity  to  inhaled  GB  in  the  higher 
species  (cat,  dog.  pig.  monkey),  the  best  possible  course  of  action  is  to  base  human  estimates  on 
available  data,  the  majority  of  which  are  derived  from  rodents.  However,  the  current  study  has 
identified  that  male  pigs  are  significantly  (p  =  0.01)  more  sensitive  to  inhaled  GB  than  female 
pigs.  While  this  conclusion  is  not,  by  itself,  enough  to  suggest  that  current  human  estimates  be 
revised,  gender  differences  in  a  species  that  more  closely  reflects  human  toxicity  estimates 
warrants  consideration,  if  not  perhaps  priority,  in  deriving  such  estimates. 

4.3  Cholinesterase  Depression  and  GB  Uptake. 

Thirty-one  of  the  32  pigs  (97%)  from  whom  the  blood  samples  were  collected 
during  exposure  showed  decreases  in  AChE  activity  to  8%  or  less  of  baseline  values  by  the  end 
of  the  exposure-duration.  Despite  the  low-levels  of  AChE  activity,  16  of  the  32  pigs  survived. 

As  early  as  1958,  Grob  and  Harvey  identified  that  red  blood  cell  cholinesterase  activity  could  be 
depressed  in  humans  to  near  zero  (with  multiple  low  dose  exposures)  without  resultant  death.27 
The  current  study  supports  the  existence  of  a  poor  correlation  between  absolute  AChE  activity 
values  and  predictability  of  lethality.  Additionally,  there  were  no  significant  differences  found 
between  the  rates  of  AChE  depression  in  animals  that  lived  versus  animals  that  died. 

Surprisingly,  rates  of  AChE  depression  were  significantly  different  between  male 
and  female  pigs,  with  the  males  showing  the  steeper  slopes.  Although  other  literature  sources 
have  not  yet  been  found  to  support  this  hypothesis,  there  are  several  possible  explanations  that  ♦ 

could  account  for  this  finding.  The  three  most  likely  are 

1)  Differences  in  the  rate  of  uptake  of  the  nerve  agent  into  the  systemic 

circulation, 

2)  Female  pigs  possessing  an  additional  buffering  capacity  within  the  circulation, 
which  acts  as  a  “sink"  to  prevent  binding  of  the  nerve  agent  to  the  cholinesterase,  and 

3)  Differences  in  the  inherent  characteristics  of  the  red  blood  cell  cholinesterase 
of  the  two  sexes. 
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In  addressing  these  possibilities,  the  evidence  provided  by  the  GB  uptake  rates  makes  differences 
in  the  breathing  rates  of  the  two  sexes  unlikely.  No  statistically  significant  differences  were 
observed  in  GB  uptake  rates  between  the  male  and  female  pigs. 

Solid  evidence  in  other  species,  especially  rodents,  supports  the  idea  of  gender 
differences  in  enzymes  with  the  potential  to  provide  a  buffering  capacity.  Female  rats  have 
higher  plasma  cholinesterase  and  carboxylesterase  activities  than  male  rats.28'30  Female  mice 
have  two- fold  higher  plasma  butylcholinesterase  activities  than  male  mice.31  In  human  adults, 
females  have  been  shown  to  have  significantly  higher  baseline  RBC  AChE  values,  but 
significantly  lower  baseline  BChE  activities  than  males.  ‘  Additionally,  there  is  evidence  that 
RBC  AChE  activity  can  be  cyclically  regulated  by  hormonal  influence  in  females.34  In  the 
present  study,  baseline  AChE  activities  in  male  and  female  pigs  show  no  statistical  differences 
(see  Section  3.5).  However,  baseline  values  of  BChE  were  significantly  higher  in  male  pigs  than 
in  female  pigs.  The  male  pigs  were  significantly  more  sensitive  to  GB.  Interestingly,  in  rats, 
females  have  higher  BChE  activity  and  they  are  also  more  sensitive  to  GB  than  their  male 
counterparts.  There  is  a  lack  of  data  available  to  determine  if  there  is  a  gender  derived  nerve 
agent  sensitivity  difference  in  humans  and  if  so  which  sex  is  more  sensitive.  It  should  be  noted, 
however,  that  like  male  pigs  and  female  rats,  men  have  higher  baseline  BChE  activities  than 
women.33 


Neither  absolute  AChE  activity  levels  nor  the  rate  of  AChE  activity  depression 
were  accurate  measurements  for  predicting  mortality.  However,  there  was  a  significant 
difference  in  the  rate  of  GB  uptake  in  animals  that  lived  versus  animals  that  died.  Not 
surprisingly,  pigs  that  died  had  a  significantly  greater  rate  of  GB  uptake  than  those  that  survived. 
Since  the  rate  of  GB  uptake  is  proportionate  to  the  concentration  of  nerve  agent,  most  likely  the 
animals  that  died  were  exposed  to  higher  concentrations.  Nonetheless,  the  data  suggest  that  the 
rate  of  nerve  agent  uptake  is  a  more  accurate  predictor  of  toxicological  endpoints  in  whole-body 
inhalation  exposures  with  nerve  agents  than  the  current  standard  of  analyzing  depression  of 
cholinesterase  activity. 


5.  CONCLUSIONS 

The  current  study  was  conducted  with  the  intent  of  estimating  lethal 
concentrations  of  the  nerve  agent  sarin  (GB)  as  a  function  of  exposure-duration  in  the  Gottingen 
minipig.  Ordinal  regression  was  used  to  fit  various  response  models  to  the  data.  The  ECT50 
<•  (severe)  and  LCT50  values  were  calculated  in  male  and  female  pigs  exposed  to  GB  vapor  for  1 0, 

60,  and  1 80  min.  The  value  of  the  toxic  load  exponent  was  essentially  independent  of  the  model 
used.  The  toxic  load  exponent  of  the  best-fit  model  (L5)  was  1 .38  (with  a  95%  confidence 
interval  of  1.24  to  1.52).  Because  this  interval  does  not  overlap  one.  Haber's  rule  is  not 
considered  an  appropriate  time  dependence  model  for  this  dataset.  Potential  curvature  in  the  data 
was  evaluated  by  inserting  a  (logT)2  term  into  the  model  and  this  term  was  found  to  be 
statistically  insignificant.  The  probit  slope  of  the  best  model  fit  (L5)  was  12.4  with  a  95%  range 
of  6.2  to  1 8.6.  The  models  were  tested  for  possible  gender  effects  and  Sex  was  found  to  be  a 
significant  term  (p  =  0.013  in  model  L5),  with  males  being  significantly  more  sensitive  than 
females.  The  ECT50  values  (severe)  were  approximately  71  to  79  %  of  the  LCT50  values  in  male 
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pigs  and  84  to  89%  of  the  LCT50  values  in  female  pigs.  The  ratios  of  severe  to  lethal 
concentrations  were  higher  in  female  minipigs  (99%  ANOVA  confidence)  indicating  that  there  is 
less  difference  between  severely  toxic  and  lethal  dosages  in  the  female  as  compared  to  male  pigs. 

Baseline  RBC  AChE  activities  showed  no  significant  differences  when  the  results 
of  male  and  female  minipigs  were  compared.  However,  female  minipigs  showed  significantly 
less  (p  =  0.022)  baseline  BChE  activity  than  male  minipigs.  Depression  of  AChE  and  BChE 
activity  (expressed  as  a  percent  of  baseline  measurements)  was  plotted  versus  exposure  dosage 
(mg.min/m3).  The  AChE  activity  was  decreased  to  8%  or  less  of  baseline  before  the  conclusion 
of  the  exposures  in  3 1  of  the  32  pigs,  with  very  little  subsequent  depression  in  AChE  activity 
after  the  conclusion  of  the  exposures.  Depression  of  BChE  was  variable  but  dropped  below  50% 
of  baseline  in  only  one  of  the  pigs.  There  were  no  significant  differences  between  the  rates  of 
depression  of  AChE  activity  in  the  pigs  that  lived  versus  the  pigs  that  died.  However,  the  rate  of 
depression  of  AChE  activity  was  significantly  higher  (p  =  0.043)  in  male  pigs  than  in  female  pigs 
at  dosages  below  7.0  mg.min/m3.  There  was  a  highly  significant  difference  (p  =  0.004)  between 
the  uptake  rates  of  GB  between  the  pigs  that  survived  for  24  hr  after  exposure  and  the  pigs  that 
eventually  died.  Not  surprisingly,  the  uptake  rate  was  higher  in  the  animals  that  eventually  died. 
The  data  suggest  that  the  nerve  agent  uptake  rate  is  a  more  accurate  predictor  of  toxicological 
endpoints  in  whole-body  inhalation  exposures  with  nerve  agents  than  the  current  standard  of 
analyzing  depression  of  cholinesterase  activity. 
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Table  1.  Durations  and  Concentrations  of  GB  Exposure  for  Male  and  Female  Pigs 


Sex 

Animal  # 

Time 

(minutes) 

Concentration 

(mg/m3) 

CT 

(mg.min/m3) 

Result  (1=moderate, 
2=severe,  3=lethal) 

Male 

41 

180 

1.00 

180.00 

3 

Male 

42 

10 

9.40 

94.00 

3 

Male 

43 

10 

6.70 

67.00 

3 

Male 

44 

10 

9.50 

95.00 

3 

Male 

45 

180 

1.34 

241.20 

3 

Male 

46 

60 

1.68 

100.80 

1 

Male 

47 

10 

7.45 

74.50 

2 

Male 

48 

60 

2.00 

120.00 

3 

Male 

49 

60 

2.50 

150.00 

3 

Male 

50 

180 

0.90 

162.00 

2 

Male 

51 

60 

1.64 

98.40 

2 

Male 

52 

180 

1.10 

198.00 

3 

Male 

53 

180 

0.80 

144.00 

1 

Male 

54 

10 

5.35 

53.50 

1 

Male 

55 

60 

1.78 

106.80 

2 

Male 

57 

60 

1.90 

114.00 

3 

Male 

58 

180 

0.99 

178.20 

2 

Male 

59 

10 

5.90 

59.00 

2 

Male 

60 

60 

1.70 

102.00 

3 

female 

61 

60 

1.82 

109.20 

1 

female 

62 

180 

0.61 

109.80 

1 

female 

63 

60 

1.49 

89.40 

3 

female 

64 

180 

0.84 

151.20 

2 

female 

65 

10 

7.95 

79.50 

1 

female 

66 

60 

1.28 

76.80 

1 

female 

67 

10 

5.10 

51.00 

1 

female 

68 

180 

1.06 

190.80 

3 

female 

69 

10 

9.74 

97.40 

1 

female 

70 

10 

6.53 

65.30 

2 

female 

71 

180 

0.91 

163.80 

1 

female 

73 

10 

10.50 

105.00 

3 

female 

74 

10 

12.78 

127.80 

3 

female 

75 

60 

1.71 

102.60 

1 

female 

76 

60 

2.52 

151.20 

2 

female 

77 

180 

1.06 

190.80 

3 

female 

78 

60 

3.08 

184.80 

3 

female 

79 

60 

2.63 

157.80 

3 

female 

80 

10 

8.00 

80.00 

3 
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Table  2.  MLE  for  Median  Effective  Concentrations  and  Dosages  (with  approximate  95% 
confidence  intervals  on  the  dosages) 


Lethality 


Males 

Females 

Exposure- 

duration 

(minutes) 

LC?o 

LCTjo 

95%  Limits 

LC5.i 

LCT50 

95%  Limits 

10 

7.25 

72.5 

55.1—95.2 

8.69 

86.9 

67.3—112.3 

60 

1.76 

105.7 

83.7—133.5 

2.12 

127.1 

98.5—163.9 

180 

1.01 

182.3 

140.6—236.3 

0.97 

174.2 

129.4—234.7 

Severe  Effects 


Males 

Females 

Exposure- 

duration 

(minutes) 

EC50 

ECTso 

95%  Limits 

EC50 

ECTS„ 

95%  Limits 

10 

5.15 

51.5 

36.9—71.9 

7.74 

77.4 

60.5—99.0 

60 

1.38 

83.0 

62.6—110.0 

1.88 

112.5 

86.6—146.0 

180 

0.74 

134.0 

97.6—182.3 

0.81 

145.9 

108.4—196.5 
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Table  3.  Ratios  of  ECT50  (severe)  Values  to  LCT50  Values 


Duration  (min)  Gender  Severe  /  lethal  Gender  Severe  /  lethal 

10 

Male 

0.711 

Female 

0.891 

60 

Male 

0.785 

Female 

0.885 

180 

Male 

0.735 

Female 

0.837 

4 


Table  4.  Probit  Slopes  and  Toxic  Load  Exponents  (n)  Obtained  from  Various  Ordinal 
Logistic  Regression  Model  Fits  without  Pig  63 


ID 

Terms  in  Model 

kc 

SE(C) 

kT 

SE(T) 

n 

SE(n) 

LI 

LogC  Time  Sex  Time*Sex 

15.7 

4.0 

— 

— 

... 

— 

L2 

LogC  Time  Sex 

12.6 

3.2 

— 

... 

— 

... 

L3 

LogC  Time 

10.2 

2.7 

— 

... 

— 

L4 

LogC  LogT  Sex  LogT*Sex 

13.5 

3.4 

9.8 

2.5 

1.38 

0.06 

L5 

LogC  LogT  Sex 

12.4 

3.1 

9.0 

2.3 

1.38 

0.07 

L6 

LogC  LogT 

9.9 

2.5 

7.2 

1.9 

1.37 

0.06 

Table  5.  Probit  Slopes  and  Toxic  Load  Exponents  (n)  Obtained  from  Various  Ordinal 
Logistic  Regression  Model  Fits  with  Pig  63 


ID 

Terms  in  Model 

kc 

SE(C) 

kT 

SE(T) 

n 

SE(n) 

LI 

LogC  Time  Sex  Time*Sex 

10.7 

2.8 

— 

... 

... 

... 

L2 

LogC  Time  Sex 

9.8 

2.6 

— 

... 

... 

— 

L3 

LogC  Time 

9.0 

2.5 

— 

— 

... 

... 

L4 

LogC  LogT  Sex  LogT*Sex 

9.9 

2.6 

7.2 

1.9 

1.37 

0.08 

L5 

LogC  LogT  Sex 

9.2 

2.4 

6.7 

1.8 

1.37 

0.08 

L6 

LogC  LogT 

8.4 

2.3 

6.2 

1.7 

1.36 

0.09 

r 
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Figure  1 .  MLE  LCT50  Estimates  as  a  Function  of  Exposure  Time 
for  Male  and  Female  Minipigs 
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Figure  2.  MLE  ECT50  (Severe)  Estimates  as  a  Function  of  Exposure  Time 
for  All  Male  and  Female  Minipigs 
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Figure  3.  Toxic  Load  Fits  (Model  L5)  of  MLE  LCT50  Estimates 
for  all  Male  and  Female  Minipigs 
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Figure  4.  Toxic  Load  Fits  (Model  L5)  of  MLE  ECT50  (Severe)  Estimates 
for  All  Male  and  Female  Minipigs 
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Table  6.  Depression  of  AChE  and  BChE  Activity 


Pig# 

Sex 

low 

%  AChE 

during 

low 

%  BChE 

during 

%  AChE 

lowest 

%  BChE 

lowest 

67 

f 

13 

100 

6.5 

100 

63 

f 

4 

80 

4 

74* 

70 

f 

3 

95 

3 

75 

64 

f 

4 

72 

4 

70 

66 

f 

2 

47 

1 

47 

68 

f 

<1 

69 

<1 

69# 

61 

f 

7 

82 

4 

78 

65 

f 

7 

95 

4 

90 

69 

f 

8 

99 

4 

94 

74 

f 

5 

74 

4 

61* 

71 

f 

4 

94 

3 

78 

76 

f 

7 

82 

4 

80 

77 

f 

4 

81 

4 

67* 

78 

f 

7 

81 

7 

70* 

80 

f 

6 

89 

6 

79* 

79 

f 

3 

85 

3 

74* 

75 

f 

8 

96 

8 

94 

46 

m 

3 

71 

1 

62 

44 

m 

3 

83 

3 

79* 

45 

m 

1 

75 

1 

75* 

41 

m 

4 

79 

4 

79* 

42 

m 

3 

84 

3 

71* 

43 

m 

3 

92 

3 

85* 

48 

m 

3 

68 

3 

68* 

53 

m 

4 

69 

4 

58 

51 

m 

4 

64 

2 

57 

54 

m 

5 

91 

5 

90 

52 

m 

4 

68 

4 

68* 

55 

m 

6 

91 

6 

91 

59 

m 

4 

88 

4 

85 

57 

m 

3 

70 

3 

70* 

60 

m 

2 

65 

2 

65* 

58a 

m 

N.A. 

N.A. 

4 

67 

Summary  table  of  AChE  and  BChE  values  at  their  lowest  during  the  sarin  exposure  and  the 
lowest  value  that  was  observed  during  or  after  the  exposure.  Pigs  that  died  during  the  time  of  the 
exposure  or  within  10  min  after  the  exposure  are  indicated  by  *.  Pigs  that  died  more  than  10  min 
after  the  exposure  concluded  are  indicated  by  #.  aBlood  samples  were  not  able  to  be  collected 
from  pig  58  until  the  pig  was  removed  from  the  chamber  after  the  exposure. 
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Table  7.  Depression  Rate  of  AChE  Activity  and  Uptake  of  GB 


RBC  cholinesterase  Depression 
(polynomial  equation  fits) 

RBC  GB  uptake 
(linear  equation  fits) 

Pig  # 

Gender 

X 

x2 

r2 

X 

r2 

80 

F 

-3.6994 

0.0307 

0.92 

0.5036 

0.9756 

67 

F 

-2.1663 

0.0071 

0.99 

0.1372 

0.9416 

70 

F 

-1.5801 

0.0016 

0.81 

0.1820 

0.9941 

65 

F 

-2.0681 

0.0105 

0.99 

0.2250 

0.9818 

69 

F 

-2.2803 

0.0135 

0.99 

0.1371 

0.9946 

74 

F 

-2.1821 

0.0122 

0.93 

0.3743 

0.9607 

77 

F 

-2.4840 

0.0163 

0.97 

0.1377 

0.9813 

71 

F 

-1.4506 

0.0063 

0.91 

0.0985 

0.9654 

68 

F 

-1.4246 

0.0054 

0.84 

0.2009 

0.9713 

64 

F 

-1.6824 

0.0080 

0.93 

0.1169 

0.9745 

63 

F 

-3.0448 

0.0259 

0.80 

0.2724 

0.9540 

66 

F 

-2.7108 

0.0197 

0.98 

0.1669 

0.9940 

61 

F 

-2.4151 

0.0167 

0.80 

0.1704 

0.8601 

79 

F 

-1.7392 

0.0081 

0.89 

0.1011 

0.9553 

76 

F 

-1.8013 

0.0092 

0.85 

0.1019 

0.8827 

78 

F 

-1.3982 

0.0054 

0.91 

0.2073 

0.8460 

75 

F 

-1.7762 

0.0090 

0.99 

0.1367 

0.9976 

55 

M 

-2.5474 

0.0168 

0.94 

0.2279 

0.9733 

57 

M 

-3.8283 

0.0372 

0.87 

0.3182 

0.9306 

46 

M 

-2.4232 

0.0165 

0.90 

0.1352 

0.9701 

48 

M 

-2.3165 

0.0142 

0.89 

0.2484 

0.9857 

51 

M 

-2.7459 

0.0199 

0.93 

0.2081 

0.9687 

60 

M 

-3.7345 

0.0366 

0.83 

0.2669 

0.9145 

44 

M 

-3.0024 

0.0234 

0.86 

0.2364 

0.9908 

42 

M 

-3.0258 

0.0239 

0.84 

0.2277 

0.9832 

43 

M 

-4.2493 

0.0465 

0.88 

0.2282 

0.9694 

54 

M 

-4.9404 

0.0606 

0.98 

0.2220 

0.9667 

59 

M 

-2.9860 

0.0202 

0.97 

0.1764 

0.9762 

52 

M 

-1.2620 

0.0047 

0.87 

0.1671 

0.9817 

53 

M 

-1.6855 

0.0840 

0.90 

0.1333 

0.9925 

45 

M 

-1.6755 

0.0080 

0.86 

0.1413 

0.9812 

41 

M 

-3.1684 

0.0280 

0.78 

0.2843 

0.9755 
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Table  8.  Pair  Wise  T-Test  Comparisons  of  Rates  of  AChE  Depression  at 
Various  Dosages.  Statistical  significance  (p<  0.05)  indicated  by  *. 


t-test  comparison 

Dosage 
(mg.min/nv ) 

P-value 

Males  vs.  females  -  all 

i 

0.012* 

Males  vs.  females  -  all 

5 

0.027* 

Males  vs.  females  -  all 

7 

0.043* 

Males  vs.  females  -  all 

8 

0.054 

Males  vs.  females  -  all 

10 

0.088 

Males  vs.  females  -  all 

15 

0.282 

Males  vs.  females  -  all 

20 

0.696 

Males  vs.  females  —  10  min 

5 

0.025* 

Males  vs.  females  -  10  min 

10 

0.026* 

Males  vs.  females  -  10  min 

15 

0.028* 

Males  vs.  females  -  10  min 

20 

0.039* 

Males  vs.  females  -  10  min 

25 

0.128 

Males  vs.  females  -  60  min 

5 

0.041* 

Males  vs.  females  -  60  min 

10 

0.037* 

Males  vs.  females  -  60  min 

20 

0.027* 

Males  vs.  females  -  60  min 

30 

0.016* 

Males  vs.  females  -  60  min 

35 

0.071 

Males  vs.  females  —  180  min 

5 

0.947 

Males  vs.  females  -  180  min 

10 

0.670 

Males  vs.  females  —  180  min 

15 

0.514 

Males  vs.  females  -  180  min 

20 

0.432 

Table  9.  Pair  Wise  T-Test  Comparisons  of  GB  Uptake  in  to  RBCs 


t-test  comparison 

P-value 

Male  vs.  female  -  all 

0.470 

Male  vs.  female  -  10  min 

0.552 

Male  vs.  female  -  60  min 

0.068 

Male  vs.  female  -  180  min 

0.340 

Lived  vs.  died  -  all 

0.004* 
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APPENDIX  A 


BINARY  AND  ORDINAL  PROBIT  MODELS  AND 
THE  METHOD  OF  MAXIMUM-LIKELIHOOD  ESTIMATION 

A1.0  INTRODUCTION 

Traditionally,  median  effective  dosages  are  determined  via  the  use  of  probit 
analysis.1 2  In  conventional  probit  analysis  of  binary  response  data,  two  parameters  are  estimated 
simultaneously  from  experimental  quantal  data  using  the  method  of  maximum  likelihood 
estimation  (MLE):3'8  the  median  effective  stress,  p;  and  the  standard  deviation  of  effective 
stresses,  a.  In  toxicology,  the  effective  stress  is  the  base  10  logarithm  of  the  effective  dosage. 
Thus,  the  base  10  logarithm  of  the  median  effective  dosage.  logioCECTso),  corresponds  to  p, 
while  the  probit  slope  equals  the  inverse  of  a.  The  probit  model  can  also  be  extended  for  use 
with  ordinal  response  data  (categorical  data  that  have  three  or  more  possible  levels  with  a  natural 
ordering  [ex.  mild,  moderate  and  severe]).1 3  For  the  ordinal  probit  model,  there  are  individual 
p's  for  each  category,  but  the  ct's  for  each  category  are  assumed  to  equal  each  other. 

The  efficiency  of  the  MLE  procedure  with  the  probit  model  is  dependent  on  the 
sample  size.  Larger  sizes  provide  unbiased  and  minimum  variance  estimates  of  both  p  and  a. 
but  this  is  not  the  case  for  small  sample  sizes.  It  has  been  shown  that  when  solving  for  both  p 
and  a  with  small  sample  sizes  that  estimates  for  p  are  unbiased  (for  all  practical  purposes),  but 
estimates  for  ct  are  biased  (with  estimates  for  a  being  too  small  on  average).3  831  Furthermore, 
the  probability  of  MLE  solution  instability  increases  as  sample  size  decreases  when  trying  to 
solve  for  both  p  and  a.  Thus,  for  small  sample  sizes,  a  more  pragmatic  approach  is  commonly 
taken  by  fixing  a  at  some  set  value  (based  on  historical  knowledge  of  the  system  under  study) 
while  solving  for  p.  This  is  an  underlying  principle  of  the  up-and-down  method  for  estimating 
median  effective  stresses/dosages.6 

The  following  are  examples  of  both  a  binary  and  ordinal  probit  model  applied  to 
the  male  pig  GB  ordinal  data  (severe  effects  and  lethality)  for  ten  minutes  exposures  from  the 
present  study.  The  probit  slope  (1  /a  or  m)  was  held  constant  throughout  the  computations. 


A2.0  BINARY  PROBIT  MODEL 

A2. 1  MLE  Algorithm  for  Binary  Probit  Model. 

For  each  trial  condition,  /,  there  is  a  likelihood,  Z,/,  of  the  observed  result 

occurring: 


(Ala) 
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loge(i, )  =  *,  [logep,]  +  («,  -  or, )  [log,  (1  -p, )] 


(Alb) 


where  p\  is  the  event  probability  for  test  condition  /,  w,  is  the  number  of  independent  trials  under 
the  i-th  condition,  and  x,  is  the  number  of  successes  in  n,.  The  likelihoods  for  all  test  conditions 
are  then  multiplied  together  to  arrive  at  the  likelihood.  L.  The  values  of  the  p\  s  that  are  most 
supported  by  the  quanta!  data  are  the  values  for  which  L  is  the  largest.  For  ease  of  calculation, 
the  natural  logarithm  of  the  likelihood  is  often  used. 

For  a  normal  distribution,/?,  is  defined  by  the  following  relations: 
p,=  f Z'  f(Z)dZ  (A2) 

J-  oo 


/(z)  = 


1 

VZt 


exp 


(A3) 
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where  Z,  is  the  standard  normal  random  variable  and  /(Z)  is  the  probability  density  function  (pdf) 
of  a  standard  normal  distribution.  In  toxicology,  the  values  of  the  individual  p\  s  are  a  function 
of  the  applied  dosages,  (C7)„  used  in  an  experiment  and  their  respective  distances  from  the 
median  effective  dosage,  ECT50.  This  is  reflected  in  the  following  definition  of  Z;: 

z,  =  =  m  {log10(CT)i  -  log10(ECT50)}  (A4) 

G 


where  s,  is  the  applied  stress  for  trial  condition  /  and  m  is  the  probit  slope  (equal  to  1/a).  The 
50%  response  level  (or  p  =  0.5)  corresponds  to  a  Z  value  of  zero.  The  MLE  estimate  of 
log(ECTso)  is  the  value  of  log(ECTso)  that  is  found  to  maximize  L  in  eq  A 1 . 

For  the  MLE  calculations,  the  first  and  second  derivatives  of  loge(T  ,)  with  respect 
to  p  are  used: 3  7-8 
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(A7) 


To  reach  convergence  at  the  value  of  log(ECT5o)  that  maximizes  logeZ.,  a  Newton- 
Raphson  (or  Newton's  Method)  algorithm  (or  similar  procedure)  can  be  used.3'7’81®  Using  the 
Newton-Raphson  method  for  the  present  system  of  equations  (eqs  A 1  b,  A6,  and  A7),  the 
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following  equation  is  used  to  determine  the  next  guess  for  log(ECT5o),  as  well  as  to  check  on 
convergence  at  the  MLE  for  log(ECTso): 


Vm*  =  - 


^loge^ 


( #  log^ 


(A8) 


where  p0  is  the  current  guess  for  p  (or  log(ECTjo)),  pnext  is  the  next  guess  for  p,  and 


x  5V  j 
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^ioge(0N 

X  > 


g2  loge^l 
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d  ioge(^,) 

dn2 


(A9) 


(A  10) 


The  first  and  second  derivatives  for  logeL  are  evaluated  at  p0.  LogeL  is  maximized  when  its  first 
derivative  with  respect  to  p  equals  zero.  Convergence  is  achieved  when  the  absolute  difference 
between  p0  and  pnext  is  less  than  a  predetermined  value. 

Thus,  the  following  algorithm  is  used  to  find  the  MLE  estimate  for  ECT50: 

(1)  Set  the  probit  slope  (m)  equal  to  some  fixed  value  for  the  duration  of  the 
algorithm. 

(2)  Make  an  initial  guess,  p0,  for  p  [or  log(ECT5o)]. 

(3)  Calculate  Z;  ,f(Zj)  and  px  for  each  test  condition  i,  corresponding  to  some  (CT), 
exposure  using  eqs  A2,  A3,  and  A4. 

(4)  Using  eq  Al,  calculate  the  individual  likelihoods,  Lj. 

(5)  Multiply  the  L,‘s  (or  add  the  loge(Z.;)'s)  together  to  estimate  the  total 
likelihood,  L  (or  logeL),  of  the  MLE  estimate. 

(6)  Calculate  the  first  and  second  derivatives  for  logeL  (evaluated  at  p0)  using 
eqs  A6,  A7,  A9,  and  A 10. 

(7)  Check  to  verify  whether  the  maximum  value  of  L  has  been  obtained.  If  not. 
go  back  to  Step  (3)  with  a  new  guess,  pnext,  for  p  [or  log(ECTso)]  using  eq  A8. 

After  the  final  log(ECT50)  estimate,  /} ,  is  obtained,  there  are  three  common  and 
general  methods  for  obtaining  approximate  confidence  limits  for  the  estimate:3  Wald  test, 
likelihood-ratio  test,  and  the  score  (or  Lagrange-multiplier)  test.  These  approximations  grow 
more  accurate  as  the  sample  size  gets  larger. 
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In  the  present  study,  the  Wald  test  was  used  to  calculate  confidence  limits.  Limits 
from  the  Wald  test  can  be  readily  obtained  from  calculations  performed  as  part  of  the  Newton- 
Raphson  algorithm  used  for  finding  the  maximum  value  for  L.  However,  the  likelihood-ratio  test 
required  additional  Newton-Raphson  algorithm  iterations. 

In  the  present  study,  the  following  equation  was  used  (based  on  the  Wald  test)  to 
calculate  the  95%  asymptotic  confidence  interval  for  p  or  the  log(ECT5o):3 


A  - 


(1.96) 


M  log,  Z,(A) 


<  p  <  p  + 


(196) 

l-</2  log,  Z.(pj 

v  rfp2 


(All) 


where  the  second  derivatives  for  logeZ.  are  evaluated  at  fi  using  eqs  A  7  and  A 1 0. 

A2.2  Example  of  Application  of  Binary  Probit  Model  with  Fixed  Probit  Slope. 

The  following  are  the  binary  data  for  the  10-min  exposures  of  the  male  pig  to  GB 
vapor.  Dosage  is  in  units  of  mg-min/m3. 


Table  Al:  Male  Pig  GB  Ordinal  Data  (10-min  exposure-duration) 


Pig 

Dosage  (CT) 

logio(CT) 

Outcome 

X, 

1 

53.5 

1.728354 

<  severe 

0 

2 

59.0 

1.770852 

severe 

0 

3 

67.0 

1.826075 

death 

1 

4 

74.5 

1.872156 

severe 

0 

5 

94.0 

1.973128 

death 

1 

6 

95.0 

1.977724 

death 

1 

For  this  example,  test  condition  /'  will  only  have  one  pig.  So.  for  eq  Al,  n  will 
equal  one  for  each  test  condition.  The  values  for  jr,  correspond  to  the  absence  or  presence  of 
lethality  in  the  exposed  pig. 

Steps  ( 1 )  and  (2):  Probit  Slope  and  Initial  Guess  for  logm(ECTso)  for  Iteration 

One 


For  Step  (1)  of  the  algorithm,  the  probit  slope  is  set  equal  to  10,  which  was  used 
as  the  step  size  for  the  up  and  down  method  employed  in  the  present  study.  For  Step  (2).  the 
initial  guess  for  the  logio(ECT50)  is  1.85304. 
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Step  (3):  Calculation  of  Z,‘s  and  p,  for  Iteration  One  Using  eqs  A2,  A3,  and  A4 


Z,  =  [1.72835)-  1.85304)]/ (1/10)= -1.24685 

=>  p\  =  0.10623 

Z2  =  [1.77085)-  1.85304)]/ (1/10)=  -0.82187 

=>  p2  =  0.20558 

Z3  =  [1.82608)-  1.85304)]/ (1/10)=  -0.26964 

=>  Pi  =  0.39372 

Z4  =  [1.87216)-  1.85304)]  / (1/10)  =  0.19117 

=>  p4  =  0.57580 

Z5  =  [1.97313)-  1.85304)]/ (1/10)=  1.20089 

=>  p5  =  0.88510 

Z6  =  [1.97772)-  1.85304)]  /  (1/10)  =  1.24685 

=>  p6  =  0.89377 

Steps  (4)  and  (5):  Calculation  of  L^s  and  L  for  Iteration  One  Using  eq  A1 

loge(Z.,)  =  (0)(loge(0.10623))  +  (l-0)(loge(l  - 

0.10623))  =  -0.11230 

loge(L2)  =  (0)  (loge(0. 20558))  +  (1  -  0)  (loge(  1  - 

0.20558))  =  -0.23014 

loge(Z.3)  =  (1)  (loge(0.39372))  +  (1  -  1)  (loge(l  - 

0.39372))  =  -0.93212 

loge(I4)  =  (0)  (loge(0.57580))  +  (1  -  0)  (loged  - 

0.57580))  =  -0.85756 

loge(I5)  =  (1)  (loge(0.885  10)) +  (1  -1)  (loged  - 

0.88510))  =  -0.12205 

loge(I6)  =  (1)  (loge(0. 89377))  +  (1  -  1 )  (loge(  1  ~ 

0.89377))  =  -0.11230 

Sum  of  the  above  loge(Z.i)'s  (or  loged))  equals  - 


2.36647. 


Step  (6):  Calculate  the  First  and  Second  Derivatives  for  load  (evaluated  at  un) 

Using  eqs  A6.  A7,  A9  and  A10 

d[loged,)]/d[no]  =  2.05162 
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0.85430. 


269.72. 


d[loge(Z,2)]/d[Ho]  =  3.58246 
d[loge(l3)]/d[Ho]  =  -9.77094 
d[loge(Z,4)]/d[M  =  9.23437 
d[loge(L5)]/d[M  =  -2.19159 
d[loge(Z,6)]/d[|ao]  =  -2.05162 
Sum  of  the  above  d [ I oge( Z- 1  )]/d [|a0]  ‘s  equals 

d2[loge(Z,,)]/d[M2  =  -29.790 


d2[l0ge(Z,2)]/d[p0]2 

=  -42.277 

d2[loge(^3)]/d[Po]2 

=  -69.125 

d2[loge(I4)]/d[p0]2 

=  -67.620 

d2[loge(£5)]/d[Po]2 

=  -31.122 

d2[l0ge(Z,6)]/d[p0]2 

=  -29.790 

Sum  of  the  above  c 

l2[loge(Z.i)]/d[p0]2  ‘s  equals- 

Step  (7):  Check  for  Convergence  on  Maximum  L  Value  and  New  Guess  for 

log(ECTsn)  for  Iteration  Two  Using  eq  A8 

After  the  first  iteration,  the  next  guess  for  p  is  found  to  equal: 

|anex,  =  1 .85304  -  (0.85430)  /  (-269.72)  =  1 .85304  -  (-0.003 1 7)  =  1 .8562 1 

Convergence  was  nearly  reached  after  the  first  iteration,  as  seen  above  with  only 
difference  of -0.003 1 7  between  p„ext  and  p0.  After  the  second  iteration,  the  difference  falls 
further  to  2.1  x  10'7.  Thus,  the  final  estimate  for  logio(ECTso)  is  1.85621  (or  (ECT50)  =  71.8  mg 
min/m3),  and  the  final  loge(Z.)  value  was  -2.365 12.  The  denominators  of  eq  A1 1  were  found  to 
equal  the  square  root  of  269.7.  With  this  value  for  the  denominators,  the  corresponding  95% 
asymptotic  confidence  interval  for  logio(ECT5o)  equals  1.73686  to  1.97556,  or  for  ECT50,  the 
interval  is  54.6  to  94.5  mg-min/m3. 
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A3.0 


ORDINAL  PROBIT  MODEL 


A3.1  MLE  Algorithm  for  Ordinal  Probit  Model 

To  model  an  ordinal  ternary  response,  eq  A1  is  modified  as  follows  for  each  trial 

condition  /: 


L,  =  Pxt  ~  P\i  ~  P2,i)" 


(A  12a) 


l09e(i, )  =  *i.,  [logePiJ  +  *2,  [logeP2,,]  +(«,  -  Xy  -  *2i, )  [loge  (1  -  pv  - p2l )]  (A  1 2b) 

where  p\,  and  p2j  are  the  event  probabilities  for  responses  of  categories  one  and  two, 
respectively,  for  test  condition  /,  is  the  number  of  independent  trials  under  the  /-th  condition, 
and  xij  and  X:j  are  the  number  of  responses  of  categories  one  and  two,  respectively,  in  /?,.  As 
before  with  the  binary  model,  the  likelihood,  L,  is  obtained  by  multiplying  together  the 
individual  likelihoods  (for  all  test  conditions). 

The  ordinal  categories  are  in  order  of  increasing  severity,  with  category  two  being 
more  severe  than  category  one.  For  this  example,  category  two  corresponds  to  lethality,  while 
category  one  corresponds  to  severe  effects.  Two  normal  distributions  are  represented  in  eq  12 
and  are  defined  as  follows  for  the  /-th  condition: 

P2,=  \Z2"f{Z)dZ  (A13) 

J- 00 


P,.I  +  P2.,=  f{Z)dZ 

J-  oo 


(A  14) 


The  values  of  the  individual  pi/s  and  pi/s  are  a  function  of  the  applied  dosages, 
(C7)„  used  in  an  experiment  and  their  respective  distances  from  their  corresponding  median 
effective  dosages:  for  category  one,  pi  or  ECT50  (severe);  and  for  category  two,  or  LCT50 
(lethality).  This  is  reflected  in  the  following  definitions  of  the  Z,' s  for  each  /-th  condition: 

Zu  =  *S'  -  m  {log10(CT),  -  log10(ECT50 )}  (A15) 

G 

Z2„  =  {S|  '  ^  =  m  {log10(CT)j  -  log10(LCT50)}  (A16) 

G 


where  p?  >  Pi  and  a  (or  m)  is  assumed  to  be  the  same  for  both  distributions.  MLE  is  now  used 
to  simultaneously  obtain  estimates  for  both  the  ECT50  (severe)  and  LCT50.  However,  some 
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modifications  are  needed  to  eqs  A6  and  A 7  to  account  for  the  first  and  second  derivatives  of 
loge(Z, ,)  with  respect  to  pi  and  p2,  respectively: 
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(A20) 


where  po, ,  equals  (1  -pu~  p2.  Also,  the  partial  derivative  of  loge(Z-  /)  with  respect  to  both  pi 

and  p2,  is  required: 


r  log  e(L,)) 
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v  dpydp2  J 
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/(Z,,)/(Z2.,) 


1,1 

Pi,! 


(A21) 


To  reach  convergence  at  the  values  for  pi  (or  ECT50  (severe))  and  p2  (or  LCT50 
(lethality))  that  maximizes  logeZ„  a  Newton- Raphson  (or  Newton's  Method)  algorithm  (or 
similar  procedure)  can  be  used.  ’ ' '  Using  the  Newton-Raphson  method  for  the  present  system 
of  equations  (eqs  A1 7  to  A2 1),  the  following  simultaneous  equations  are  used  to  determine  the 
next  guess  for  the  median  effective  dosages,  as  well  as  to  check  on  the  convergence  of  the 
solution: 


<?loge(Z,)l 

(  a2loge(Z,)l 

A/y,  + 

f  d2  loge(Z,)>| 

J 

l  J 

v  dp,dp2  J 

3loge(Z,)l 

d2  loge(L)^ 

A/y,  + 

'  d2  loge(Z.)^| 

dp2  ) 

,  dPydp2  j 

l  ) 

(A22) 


(A23) 


where  Api  =  (pnext  -  p0)i  and  Ap2  =  (gnext  -  Po)2  for  median  effective  dosages  1  (severe)  and  2 
(lethal),  respectively,  and 
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(A24) 
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(A26) 


where.y  equals  1  and  2  for  severe  and  lethality,  respectively.  The  above  derivatives  for  logeT  are 
evaluated  at  pi.0  and  p2.0.  LogeT  is  maximized  when  its  first  derivatives  with  respect  to  the  two 
p's  equal  zero.  Convergence  is  achieved  when  the  absolute  difference  between  p0  and  p„ext  is 
less  than  a  predetermined  value.  After  convergence  is  reached,  the  asymptotic  variance- 
covariance  matrix  of  the  maximum  likelihood  estimate  can  be  calculated  by  taking  the  inverse  of 
the  matrix  of  second  derivatives  of  LogeZ.:  ' ' 


t(a1,a2)  = 


js2logeU)' 

\  d2  i°ge  (T)l 

l  j 

,  5//,5/y2  j 

fa2loge(^ 

r52loge(L)^ 

\dp,dp2  j 

l  ^a22  J 

(A27) 


Thus,  the  following  algorithm  is  used  to  find  the  MLE  estimate  for  ECT50: 

(1)  Set  the  probit  slope  ( m )  equal  to  some  fixed  value  for  the  duration  of  the 
algorithm. 

(2)  Make  initial  guesses  for  log(ECTso)  and  log(LCT50):  pi,0  and  p2,o< 
respectively. 

(3)  Calculate  ZtJ  and  p  / j  for  each  test  condition  i  and  mean  dosage  /, 

corresponding  to  some  (CT)j  using  eqs  A3  and  A 13  to  A 16. 

(4)  Using  eq  A 12,  calculate  the  individual  likelihoods,  L,. 

(5)  Multiple  the  Lf  s  [or  add  the  loge(T,)'s]  together  to  estimate  the  total 
likelihood,  L  (or  logeZ.),  of  the  MLE  estimate. 
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(6)  Calculate  the  first  and  second  derivatives  for  logeZ,  (evaluated  at  p0)  for  each 
mean  dosage  j,  using  eqs  A17  to  A20.  Also,  calculate  the  derivative  for  loge7 
with  respect  to  both  pi,0  and  p2.o  using  eq  A2 1 . 

(7)  Check  to  verify  whether  the  maximum  value  of  L  has  been  obtained.  If  not. 
go  back  to  Step  (3)  with  new  guesses.  pi.next  and  p2jiext,  for  pi  and  |.i2, 
respectively,  by  solving  eqs  A22  and  A23,  simultaneously. 

In  some  instances  with  Step  (7).  poor  initial  guesses  for  pi  and  p2  may  produce 
the  situation  where  ja i  .next  >  P2.next.  a  violation  of  a  key  boundary  condition.  A  simple  resolution 
for  this  problem  is  to  artificially  reduce  (for  this  iteration)  the  values  of  Api  and  Ap2  by  some  set 
factor. 


After  the  final  log(ECT5o)  and  log(LCT5o)  estimates,  //,  and  fj2  .  are  obtained, 

there  are  three  common  and  general  methods  for  obtaining  approximate  confidence  limits  for 
these  estimates:3  Wald  test,  likelihood-ratio  test,  and  the  score  (or  Lagrange-multiplier)  test.  As 
with  the  binary  probit  model  (see  Section  A2.0),  these  approximations  grow  more  accurate  as  the 
sample  size  gets  larger. 

In  the  present  study,  the  Wald  test  was  used  to  calculate  confidence  limits.  Limits 
from  the  Wald  test  can  be  readily  obtained  from  calculations  performed  as  part  of  the  Newton- 
Raphson  algorithm  used  for  finding  the  maximum  value  for  L.  However,  the  likelihood-ratio  test 
required  additional  Newton-Raphson  algorithm  iterations. 

In  the  present  study,  the  following  equation  was  used  (based  on  the  Wald  test)  to 
calculate  the  95%  asymptotic  confidence  interval  for  each  m  (with  j  =  1  for  log(ECT50)  and  2  for 
log(LCT50)):3'7'8 


(196) 

Vvar(ft,) 


< 


< 


A  j 


(1-96) 

>/var(ft;) 


(A28) 


where  var (juf)  is  the  variance  for  p7.  The  values  for  the  variances  and  covariance  are  calculated 
using  eq  A27. 

A3. 2  Example  of  Application  of  Ordinal  Probit  Model  with  Fixed  Probit  Slope 

The  following  are  the  quantal  data  for  the  ten-minute  exposures  of  the  male  pig  to 
GB  vapor.  Dosage  is  in  units  of  mg-min/m3. 
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Table  A2:  Male  Pig  GB  Ordinal  Data  (10-min  exposure-duration) 


Pig 

Dosage  (CT) 

logio(CT) 

Outcome 

Xl 

*2 

1 

53.5 

1.728354 

<  severe 

0 

0 

2 

59.0 

1.770852 

severe 

1 

0 

3 

67.0 

1.826075 

death 

0 

1 

4 

74.5 

1.872156 

severe 

1 

0 

5 

94.0 

1.973128 

death 

0 

1 

6 

95.0 

1.977724 

death 

0 

1 

For  this  example,  test  condition  i  will  only  have  one  pig.  So.  n  will  equal  one  for 
each  test  condition  for  eq  A 12.  The  values  for  X/  and  x2  correspond  to  the  absence  or  presence  of 
a  maximum  effect:  x/  equals  one  if  the  maximum  effect  observed  was  severe  effects  (and  equals 
zero  otherwise);  and  x2  equals  one  if  the  maximum  effect  was  lethality  (and  equals  zero 
otherwise). 


Steps  ( 1 )  and  (2):  Probit  Slope  and  Initial  Guesses  for  login(ECTsn)  and 

login(LCTsn)  for  Iteration  One 

For  Step  (1)  of  the  algorithm,  the  probit  slope  is  set  equal  to  10,  which  was  used 
as  the  step  size  for  the  up  and  down  method  employed  in  the  present  study.  For  Step  (2),  the 
initial  guesses  for  logio(ECT5o)  and  logio(LCTso)  are  1.68  and  1.83,  respectively. 

Steps  (3),  (4)  and  (5):  Calculation  of  Z,  , ,  f(Zjj)  and  p  ,  ,  Using  eqs  A3  and  A13  to 

A16  and  of  L£s  and  L  Using  eq  A12  for  Iteration  One 

Table  A3  shows  the  calculated  values  for  each  pig.  The  sum  of  the  individual 
logeZ,,  equals  -3.8212. 


Table  A3:  Values  of  Z,t/ ,  f(ZKj)  and  p  tj  for  First  Iteration 


Pig 

z, 

z2 

f(Zi) 

f(Z2) 

pi 

P2 

Loge  Lj 

1 

0.48354 

-1.01646 

0.35493 

0.23799 

0.53094 

0.15470 

-1.15720 

2 

0.90852 

-0.59148 

0.26404 

0.33492 

0.54110 

0.27710 

-0.61420 

3 

1.46075 

-0.03925 

0.13727 

0.39864 

0.44361 

0.48434 

-0.72500 

4 

1.92156 

0.42156 

0.06297 

0.36502 

0.30934 

0.66333 

-1.17330 

5 

2.93128 

1.43128 

0.00543 

0.14324 

0.07449 

0.92382 

-0.07920 

6 

2.97724 

1.47724 

0.00474 

0.13398 

0.06835 

0.93019 

-0.07240 
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Step  (6):  Calculate  the  Various  Derivatives  for  logPZ,  (evaluated  at  ui  n  and  u?_n) 
for  Each  Mean  Dosage  /,  Using  eqs  A17  to  A21. 

Table  A4  shows  the  calculated  values  for  each  pig.  Pigs  showing  severe  effects 
make  contributions  to  all  five  derivatives,  whereas  those  pigs  having  less  than  severe  effects  or 
died  only  contribute  to  two  of  the  five  derivatives. 


Table  A4:  Values  of  logPZ,  Derivatives  for  First  Iteration 


Pig 

(d\og  e(Z,n 

l  J 

f32loge(L,)^ 

l  J 

l  d//2  J 

p2ioge(i,)^ 

l  S//22  J 

( <32  loge(L,)Nj 

[  dp,  dp2  J 

1 

11.2906 

-72.8829 

0 

0 

0 

2 

-4.8798 

-68.1456 

6.1896 

-74.9220 

30.2039 

3 

0 

0 

-8.2304 

-64.5089 

0 

4 

-2.0355 

-43.2572 

11.8000 

-89.4952 

24.0191 

5 

0 

0 

-1.5505 

-24.5966 

0 

6 

0 

0 

-1.4404 

-23.3522 

0 

Total 

4.3753 

-184.2857 

6.7683 

-276.8749 

54.2231 

Step  (7):  Check  for  Convergence  on  Maximum  L  Value  and  New  Guesses  for 

log(ECTsn)  and  log(LCTsn)  for  Iteration  Two  Using  eqs  A22  and  A23. 

Solving  eqs  A22  and  A23  simultaneously  produces  values  of  0.03283  and 
0.03087  for  Api  and  Api,  respectively.  So,  the  next  guesses  for  pi  and  p2  are  1 .71283  and 
1.86087,  respectively.  After  the  third  iteration,  the  absolute  values  for  both  Api  and  Ap;  are  less 
than  1  x  10°.  At  which  point,  the  final  values  for  pi  and  pi  are  1.71203  and  1.86015, 
respectively.  Thus,  the  final  estimates  for  ECT50  and  LCT50  are  51.5  and  72.5  mg-min/m3, 
respectively.  The  final  loge(L)  value  was  -3.6478.  The  final  variance-covariance  matrix  equals: 

0.005449  0.001100' 

0.001100  0.003642 

Using  the  above  information  with  eq  A28,  the  corresponding  95%  asymptotic 
confidence  interval  for  logi0(ECT50)  equals  1.56735  to  1.85671,  or  for  ECT50,  the  interval  is  36.9 
to71.9  mg-min/m3.  For  lethality,  the  confidence  interval  for  logio(LCT5o)  equals  1.74186  to 
1.97845,  or  for  LCT50,  the  interval  is  55.2  to  95.2  mg-min/m3. 


f(ai,A2)  = 


A4.0  OTHER  USEFUL  RELATIONSHIPS 

For  simplicity,  many  of  the  derivative  equations  are  presented  without  benefit  of 
showing  intermediate  steps.  The  following  is  a  listing  of  relationships  that  were  useful  in 
arriving  at  the  final  relations. 
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A4.1 


Binary'  Probit  Model 


A4.2 


-/(*,) 

d/J  a 

dv 


f{Z)Z 


(A29) 


(A30) 


dZ^ 

dp 


Ordinal  Probit  Model 


<K,  =  -/(Z,„) 

dp,  O 

dP2.i 


dp, 

dPi,, 

dp2 


=  0 


d/>2,, 

^//2  , 


f(Z 2„) 


dp2t,  =  ~/(Z2.,) 

d//2  cr 

g/(zi)  r  r/(zozi 

dp. 


S/(Z2)  _  <y(z,) 


az, 

dp. 


az2 

dp  2 


dp, 

iT 

yO-J 


dp  2 


=  0 


<T(Z2) 

dp2 


/(^ 


(A3 1) 

(A32) 

(A33) 

(A34) 

(A35) 

(A36) 

(A37) 

(A38) 

(A39) 


az,  az, 

— -  =  — -  =  o 

dp  2  dp, 


(A40) 
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A5.0 
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APPENDIX  B 


MINITAB  PRINTOUTS  (ALL  PIGS) 


Probit  Analysis:  Death  versus  CT,  Group 

Distribution:  Lognormal  base  10 
Response  Information 


Var iable 

Value 

Count 

Death 

1 

18  (Event 

) 

0 

20 

Total 

38 

Factor  Ii 

nformat ion 

Factor 

Levels  Values 

Group 

6  M10 

M60  Ml 80  F10 

F60 

FI  80 

Estimation  Method:  Maximum  Likelihood 

Regression  Table 

Standard 

Variable 

Coef 

Error 

z 

P 

Constant 

20.074 

6.099 

-3.29 

0.001 

CT 

10.816 

3.275 

3.30 

0.001 

Group 

M60 

-1.825 

1  .005 

1.82 

0.069 

Ml  80 

-4.335 

1.566 

-2.77 

0.006 

F10 

- 1 .0542 

0.9125 

1  .  16 

0.248 

F60 

-2.572 

1  .  105 

-2.33 

0.020 

FI  80 

-4 . 163 

1.545 

-2.69 

0.007 

Natura 1 

Response 

0.000 

Test  for 

equa 1  s 1  opes 

:  Chi-Square 

=  5.7971 .  1 

Log - L i ke 1 i hood  =  -17. 

891 

Multiple 

degree  of  freedom  test 

Term 

Ch i -Square 

DF  P 

Group 

8.823 

5  0.116 

DF  =  5.  P- Value  =  0.326 


Goodness -of -Fi t  Tests 


Method 

Ch i -  Square 

DF 

P 

Pearson 

39.051 

30 

0.  125 

Dev i ance 

35.783 

30 

0.215 

Hosmer  Lemeshow 

5.002 

8 

0.757 

Table  of  Observed  and  Expected  Frequencies: 

(See  Hosmer  Lemeshow  Test  for  the  Pearson  Chi-Square  Statistic) 


Value 

1 

Obs 

Exp 

0 


Group 
5  6 


0 

0.0 


0.4 


0 

0.9 


.2 


.5 


3 

3.2 


3 

2.2 


3 

3.5 


10  Total 

4  18 

3.8 
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Obs 

Exp 


3 

3.0 


3 

3.6 


4 

3.  1 


3 

2.8 


2 

2.4 


2 

1  .5 


2 

1  .8 


0  1  0  20 

0.8  0.5  0.2 


Total  3  4 


4  4  4 


3  5  3 


4  4  38 


Group  =  M 1 0 

Table  of  Percentiles 

Standard 

Percent  Percentile  Error 

50  71.7733  9.4814 


95.0%  Fiducial  Cl 
Lower  Upper 

52.2478  99.4268 


Group  =  M60 

Table  of  Percentiles 

Standard 

Percent  Percentile  Error 

50  105.8463  11.7477 


95.0%  Fiducial  Cl 
Lower  Upper 

79.7964  137.0774 


Group  =  M 1 80 
Table  of  Percentiles 

Standard 

Percent  Percentile  Error 

50  180.6076  22.4920 


95.0%  Fiducial  Cl 
Lower  Upper 

133.8512  245.4973 


Group  =  F10 

Table  of  Percentiles 

Standard 

Percent  Percentile  Error 

50  89.8318  11.3579 


95.0%  Fiducial  Cl 
Lower  Upper 

66.5600  123.2342 


Group  =  F60 

Table  of  Percentiles 

Standard 

Percent  Percentile  Error 

50  124.0963  15.1873 


95.0%  Fiducial  Cl 
Lower  Upper 

94.6559  172.2615 


Group  =  F180 


Table  of  Percent 

i  les 

Standard 

95 . 0% 

Fiducial  Cl 

Percent  Percentile 

Error 

Lower 

Upper 

50  174. 

1303  25 

.0783  1 

23.9176 

249.9727 

Table  of  Relative  Potency 

Factor:  Group 

Relat i ve 

95 . 0% 

Fiduc ia 

1  Cl 

Comparison 

Potency 

Lower 

Upper 

M 10  VS  M60 

1.4747 

0.9528 

2 

.2100 

M10  VS  Ml 80 

2.5164 

1 .6165 

3 

.9131 

M 10  VS  F10 

1.2516 

0.8051 

1 

.9612  not  d i f ferent 

M10  VS  F60 

1 . 7290 

1 . 1425 

2 

.7473 
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M 10 

VS 

FI  80 

2.4261 

1 .5174 

3.9296 

M60 

VS 

M 1 80 

1.7063 

1 . 1541 

2.6031 

M60 

VS 

F10 

0.8487 

0.5743 

1 . 3056 

M60 

VS 

F60 

1  . 1724 

0.8136 

1 .8322  not 

d i f  ferent 

M60 

VS 

FI  80 

1.6451 

1.0796 

2.6231 

M 1 80 

VS 

F10 

0.4974 

0.3243 

0.7696 

M 1 80 

VS 

F60 

0.6871 

0.4603 

1 .0780 

Ml  80 

VS 

FI  80 

0.9641 

0.6108 

1 . 5434  not 

d i f  ferent 

F10 

VS 

F60 

1.3814 

0.9190 

2. 1630 

F10 

VS 

FI  80 

1.9384 

1.2190 

3.0979 

F60 

VS 

FI  80 

1.4032 

0.8707 

2. 1819 

MTB  >  code  (1)0  (2  3)1  'Score'  'Severe' 
MTB  >  Probit  'Severe'  =  ' CT'  'Group': 
SUBC>  Factors  Group' ; 

SUBC>  Lgnten; 

SUBC>  Brief  3; 

SUBC>  Con  f i dence  95.0. 


Probit  Analysis:  Severe  versus  CT,  Group 

Distribution:  Lognormal  base  10 

Response  Information 

Variable  Value  Count 

Severe  1  27  (Event) 

0  1 1 

Total  38 

Factor  Information 


Factor  Levels  Values 

Group  6  M10  M60  M180  F10  F60  F180 

Es  t i ma t i on  Me  t  hod :  Max i mum  L i ke 1 i hood 


Regression  Table 


Standard 

Var iable 

Coef 

Error 

Constant 

- 16.783 

6.417 

CT 

9.793 

3.568 

Group 

M60 

-2.038 

1 .285 

Ml  80 

-4.095 

1 .852 

F10 

-1.776 

1 .024 

F60 

3.119 

1 .279 

FI  80 

-4.483 

1 .789 

Natural 

Response 

0.000 

Test  for 

equal  slopes: 

Ch i -Square 

Log -Like! 

l i hood  =  15.1 

195 

Multiple 

degree  of  freedom  test 

Term 

Ch i -  Square 

DF  P 

Group 

7.790 

5  0.168 

Goodness - of  - Fi t  Tests 


Method  Chi-Square 


Z  P 
-2.62  0.009 
2.74  0.006 

-1.59  0. 1 13 
-2.21  0.027 
1.73  0.083 
2.44  0.015 
-2.51  0.012 


=  3.7545.  DF  =  5.  P-Value  = 


DF  P 


.585 
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Pearson 
Deviance 
Hosmer - Lemeshow 


27.976 

30.389 

5.001 


30  0.572 

30  0.446 

8  0.757 


Table  of  Observed  and  Expected  Frequencies: 

(See  Hosmer  Lemeshow  Test  for  the  Pearson  Chi  Square  Statistic) 

Group 


Value 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Total 

i 

Obs 

0 

2 

2 

2 

2 

3 

4 

4 

4 

4 

27 

Exp 

0.2 

1  .4 

2.  1 

2.7 

3.2 

2.6 

3. 

.5 

3. 

7 

3.9 

4 

.0 

0 

Obs 

3 

2 

2 

2 

2 

0 

0 

0 

0 

0 

1  1 

Exp 

2.8 

2.6 

1.9 

1  .3 

0.8 

0.4 

0. 

.5 

0. 

3 

0.  1 

0. 

0 

Total 

3 

4 

4 

4 

4 

3 

4 

4 

4 

4 

38 

Group  =  M10 

Table  of  Percentiles 

Percent  Percent i le 
50  51.7458 


Standard  95.0%  Fiducial  Cl 

Error  Lower  Upper 

9.6069  25.1350  75.2237 


Group  =  M60 

Table  of  Percentiles 

Percent  Percentile 
50  83.5637 


Standard  95.0%  Fiducial  Cl 

Error  Lower  Upper 

14.1420  38.7730  111.9375 


Group  =  Ml  80 

Table  of  Percentiles 

Percent  Percentile 
50  135.5373 


Standard  95.0%  Fiducial  Cl 

Error  Lower  Upper 

24.5361  63.7336  190.5009 


Group  =  F10 

Table  of  Percentiles 

Percent  Percent i 1 e 
50  78.5629 


Standard  95.0%  Fiducial  Cl 

Error  Lower  Upper 

10.0979  52.7437  108.7298 


Group  =  F60 

Table  of  Percentiles 

Percent  Percentile 
50  107.7515 


Standard  95.0%  Fiducial  Cl 

Error  Lower  Upper 

15.2617  73.9244  163.5522 


Group  =  FI  80 
Table  of  Percentiles 


Standard 


95.0%  Fiducial  Cl 
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Error 

23.2337 


Lower 
88. 1604 


Upper 

214.5667 


Percent  Percent i  le 
50  148.4678 


Table  of  Relative  Potency 


Factor:  Group 


Relative 

95 . 0% 

Fiducial  Cl 

Compar i son 

Potency 

Lower 

Upper 

M 10 

VS  M60 

1 .6149 

0.7930 

2.8947 

M 10 

VS  Ml 80 

2.6193 

1 .2872 

4 . 9888 

M10 

VS  F10 

1.5182 

0.9208 

3.2940 

M10 

VS  F60 

2.0823 

1 .2722 

5.0262 

M10 

VS  FI 80 

2.8692 

1 .6220 

6. 1679 

M60 

VS  M 1 80 

1.6220 

0.8974 

3.1172 

M60 

VS  F10 

0.9402 

0.6218 

2. 1251 

M60 

VS  F60 

1 . 2895 

0.8486 

3.2828 

M60 

VS  FI 80 

1 .7767 

1 . 1027 

3.9527 

Ml  80 

VS  F10 

0.5796 

0.3649 

1.2943 

Ml  80 

VS  F60 

0.7950 

0.5014 

1 .9863 

Ml  80 

VS  FI 80 

1 .0954 

0.6452 

2.4146 

F10 

VS  F60 

1.3715 

0.8463 

2.4912 

F10 

VS  FI 80 

1 . 8898 

1 .0337 

3. 1909 

F60 

VS  FI 80 

1 .3779 

0.6856 

2.2820 

di f ferent 


d i f  ferent 


d i f  ferent 


MTB  >  OLogistic  ’Score'  =  logC  Time  Gender  Time*Gender; 
SUBC>  Factors  ’Time’  ’Gender’: 

SUBC>  Normit; 

SUBC>  Brief  3. 


Model  LI 

Ordinal  Logistic  Regression:  Score  versus  logC,  Time,  Gender 

Link  Function:  Normit 

Response  Information 

Variable  Value  Count 

Score  1  1 1 

2  9 

3  18 

Total  38 

Factor  Information 

Factor  Levels  Values 
Time  3  10  60  180 

Gender  2  female  male 


Logistic  Regression  Table 


Predictor 

Coef 

SE  Coef 

Z  P 

Const ( 1 ) 

9.344 

2.574 

3.63  0.000 

Const (2) 

10.268 

2.646 

3.88  0.000 

logC 

-10.698 

2.765 

-3.87  0.000 

Time 

60 

-6.624 

1.891 

-3.50  0.000 

180 

-10.363 

2.776 

-3.73  0.000 

Gender 

male 

-1 .3012 

0.7854 

-1.66  0.098 
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Time*Gender 
60*male 
180*ma 1 e 


0.248 
1 .287 


1  .031 
1  .  104 


0.24  0.810 
1  .  17  0.243 


Tests  for  terms  with  more 

Term  Chi-Square 
Time  13.936 
Time*Gender  1 . 547 


than  1  degree  of  freedom 

DF  P 

2  0.001 
2  0.461 


Log- l i kel ihood  =  29.365 

Test  that  all  slopes  are  zero:  G  =  21.369.  DF  =  6 .  P- Value  =  0.002 


Goodness  - o f - F i t  Tests 


Me  t  hod 

Chi -Square 

DF 

P 

Pearson 

77.645 

66  0 

.  155 

Dev i a nee 

58.731 

66  0 

.725 

Measures  of 

Assoc iat ion 

(Between  the  Response  Variable  and  1 

Predicted  Probabilities) 

Pa  i  rs 

Number 

Percent 

Summary  Measures 

Concordant 

390 

85 . 0% 

Somers'  D 

0.70 

Discordant 

67 

1 4 . 6% 

Goodman -Kruskal  Gamma 

0.71 

Ties 

2 

0 . 4% 

Kendal  1 ’ s  Tau  -  a 

0.46 

Tota  1 

459 

100.0% 

MTB  >  OLogistic  'Score'  =  logC  Time  Gender 
SUBC>  Factors  'Time'  ’Gender'; 

SUBC>  Normit; 

SUBC>  Brief  3. 


Model  L2 


Ordinal  Logistic  Regression:  Score  versus  logC,  Time,  Gender 


Link  Function:  Normit 

Response  Information 

Variable  Value  Count 

Score  1  1  1 

2  9 

3  18 

Total  38 


Factor  Information 

Factor  Levels  Values 
Time  3  10  60  180 

Gender  2  female  male 


Logistic  Regression  Table 


Predictor  Coef 
Const  (1)  8.353 
Const  (2)  9.250 
logC  -9.825 
Time 


SE  Coef  Z  P 

2.321  3.60  0.000 

2.389  3.87  0.000 

2.591  -3.79  0.000 


APPENDIX  B 


52 


-5.992 

-8.904 


1.646 

2.394 


-3.64  0.000 
-3.72  0.000 


i 


♦ 


60 
180 

Gender 

male  -0.7947  0.4274  -1.860.063 


Tests  for  terms  with  more  than  1  degree  of  freedom 


Term  Chi-Square  DF  P 

Time  13.880  2  0.001 


Log  1 i ke 1 i hood  =  -30.139 

Test  that  all  slopes  are  zero:  G  =  19.821.  DF  =  4 .  P-Value  =  0.001 
Goodness -of -Fi t  Tests 


Method  Chi-Square 
Pearson  71.098 
Deviance  60.278 


DF  P 

68  0.375 

68  0.736 


Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Pa  i  rs 

Number 

Percent 

Summary  Measures 

Concordant 

377 

82. 1% 

Somers’  D 

0.64 

Discordant 

82 

17.9% 

Goodman- Kruska 1  Gamma 

0.64 

Ties 

0 

0 . 0% 

Kendal  1 ’ s  Tau-a 

0.42 

Total 

459 

100.0% 

MTB  >  OLogistic  ’Score’  =  logC  Time 
SUBC>  Factors  'Time'  : 

SUBC>  Normit; 

SUBC>  Brief  3. 


Model  L3 


Ordinal  Logistic  Regression:  Score  versus  logC,  Time 


Link  Function:  Normit 

Response 

Informat  ion 

Variable 

Value  Count 

Score 

1 

11 

2 

9 

3 

18 

Total 

38 

Factor  Information 

Factor 

Levels  Values 

T  i  me 

3  10  60 

180 

Log i s  t i c 

Regression  Tabl 

le 

Predictor 

Coef 

SE  Coef 

Z 

P 

Const ( 1 ) 

7.292 

2.156 

3.38 

0.001 

Const (2) 

8.  125 

2.212 

3.67 

0.000 

logC 

-9.016 

2.463 

-3.66 

0.000 

Time 

60 

-5.540 

1.587 

-3.49 

0.000 

180 

-8.220 

2.278 

-3.61 

0.000 
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Tests  for  terms  with  more  than  1  degree  of  freedom 


Term  Chi -Square  DF  P 

Time  13.020  2  0.001 


Log  1  i ke 1 i hood  =  -31.919 

Test  that  all  slopes  are  zero:  G  =  16.261.  DF  =  3.  P  Value  =  0.001 
Goodness  -  of  - Fi t  Tests 


Method 
Pearson 
Dev i ance 


Chi -Square 
75.833 
63.839 


DF  P 

69  0.268 

69  0.653 


Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Pairs 

Number 

Percent 

Summary  Measures 

Concordant 

361 

78 . 6% 

Somers'  D 

0.58 

Discordant 

97 

21 . 1% 

Goodman - Kruska 1  Gamma 

0.58 

Ties 

1 

0.2% 

Kendal  1 's  Tau-a 

0.38 

Total 

459 

100.0% 

MTB  >  Name  m2  =  ' XPWX2 ’ 

MTB  >  OLogistic  'Score'  =  logC  logT  Sex  Sex*logT; 
SUBC>  Normit; 

SUBC>  XPWX i nverse  ' XPWX2 ' ; 

SUBC>  Brief  3. 


Model  L4 

Ordinal  Logistic  Regression:  Score  versus  logC,  logT,  Sex 


Link  Function:  Normit 

Response  Information 

Variable  Value  Count 

Score  1  1 1 

2  9 

3  18 

Total  38 


Logistic  Regression  Table 


Predictor 

Coef 

SE  Coef 

Cons  t  ( 1 ) 

15.063 

4  .  144 

Const (2) 

15.958 

4.212 

logC 

-9.865 

2.590 

logT 

-7.207 

1.924 

Sex 

-  1 . 1428 

0.7731 

logT* Sex 

0.4457 

0.4342 

Log- 1  ike  1 i hood  =  -30.147 
Test  that  all  slopes  are  zero:  G  = 

Goodness  - of  - Fi t  Tests 


Z  P 
3.63  0.000 
3.79  0.000 
-3.81  0.000 
-3.74  0.000 
1.48  0. 139 
1.03  0.305 


19.805.  DF  =  4.  P-Value  =  0.001 


Method  Chi-Square  DF  P 

Pearson  78.211  68  0.186 
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Dev i ance 


60.294 


68  0.736 


Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Pairs 

Number 

Percent 

Summary 

Measures 

Concordant 

382 

83 . 2% 

Somers ' 

D 

0.67 

Discordant 

76 

16.6% 

Goodman 

-Kruskal 

Gamma 

0.67 

Ties 

1 

0 . 2% 

Kenda 1 1 

‘s  Tau-a 

0.44 

Total 

459 

100.0% 

1 

MTB  >  print 

m2 

♦ 

Data  Display 

Matrix  XPWX2 

17. 1740 

17.4202 

_ 

10.5623 

-7.9504 

1 .0833 

0.5123 

17.4202 

17.7383 

- 

10.7412 

-8.0818 

1 . 1036 

0.5202 

10.5623 

10.7412 

6.7076 

4.8670 

0.7074 

-0.3337 

-7.9504 

8.0818 

4.8670 

3.7035 

0.4997 

-0.2371 

- 1 .0833 

-  1 . 1036 

0.7074 

0.4997 

0.5977 

-0.3224 

0.5123 

0.5202 

-0.3337 

-0.2371 

0.3224 

0.1885 

MTB  >  Name  m3  =  ' XPWX3 ' 

MTB  >  OLogistic  'Score'  =  logC  logT  Sex 
SUBC>  Normit; 

SUBC>  XPWX inverse  'XPWX3' ; 

SUBC>  Brief  3. 


Model  L5 


Ordinal  Logistic  Regression:  Score  versus  logC,  logT,  Sex 


Link  Function:  Normit 


Response 

Informat  ion 

Variable 

Value 

Count 

Score 

1 

1 1 

2 

9 

3 

18 

Total 

38 

Log i st ic 

Regression  Table 

Predictor 

Coef 

SE  Coef 

Const (1) 

14.004 

3.926 

Const (2) 

14.877 

3.990 

logC 

-9. 181 

2.448 

logT 

-6.711 

1.824 

Sex 

-0.3901 

0.2126 

Log- 1  ike  1 i hood  =  -30.692 

Test  that  all  slopes  are  zero:  G 

Goodness -of -Fi t  Tests 


Z  P 
3.57  0.000 
3.73  0.000 
-3.75  0.000 
3.68  0.000 
-1.83  0.067 


18.716.  DF  =  3.  P  Value  =  0.000 


Method  Chi-Square  DF  P 

Pearson  79.671  69  0.178 
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Deviance 


61.384 


69  0.731 


Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Pa  i  rs 

Number 

Percent 

Summary  Measures 

Concordant 

374 

8 1  . 5% 

Somers’  D 

0. 

63 

Discordant 

84 

18.3% 

Goodman  Kruska 1  Gamma 

0. 

63 

Ties 

1 

0 . 2% 

Kenda 1 1 ' s  Tau-a 

0. 

41 

Total 

459 

100.0% 

MTB  >  print  m3 

Data  Display 


Matrix  XPWX3 


5.4117 

15.6317 

-9.4455 

-7.1341 

-0. 1913 

5.6317 

15.9208 

-9.6069 

-7.2526 

-0. 1981 

9 . 4455 

-9.6069 

5.9921 

4.3492 

0. 1239 

7.1341 

-7.2526 

4.3492 

3.3254 

0.0879 

0. 1913 

-0. 1981 

0.1239 

0.0879 

0.0452 

MTB  >  Name  m4  =  ' XPWX4 ' 

MTB  >  OLogistic  'Score'  =  logC  logT 
SUBC>  Normit; 

SUBC>  XPWX inverse  ' XPWX4  '  ; 

SUBC>  Brief  3. 


Model  L6 


Ordinal  Logistic  Regression:  Score  versus  logC,  logT 

Link  Function:  Normit 


Response 

Informat  ion 

Var i ab le 

Va  1  ue 

Count 

Score 

1 

1 1 

2 

9 

3 

18 

Total 

38 

Log i s  t i c 

Regression  Table 

Predictor 

Coef 

SE  Coef 

Cons  t ( 1 ) 

12.838 

3.690 

Const (2) 

13.651 

3.747 

l°gC 

-8.402 

2.303 

logT 

-6. 180 

1 .719 

Log- 1 i ke 1 

i hood  =  -  32 

.428 

Test  that 

all  si  opes 

are  zero:  < 

Goodness- of -Fi t  Tests 


Method 

Chi -Square 

DF 

Pearson 

78.460 

70 

Dev i ance 

64.856 

70 

Z  P 
3.48  0.001 
3.64  0.000 
-3.65  0.000 
-3.59  0.000 


=  15.244.  DF  =  2.  P-Value  =  0.000 


228 

651 
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Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Pai  rs 

Number 

Percent 

Summary  Measures 

Concordant 

356 

77.6% 

Somers '  D 

0.56 

Discordant 

101 

22 . 0% 

Goodman -Kruskal  Gamma 

0.56 

Ties 

2 

0.4% 

Kendal  1 ' s  Tau  -  a 

0.36 

Total 

459 

100.0% 

MTB  >  print  m4 

Data  Display 

Matrix  XPWX4 


13.6191 

13.7995 

-8.3393 

-6.3225 

13.7995 

14.0401 

-8.4723 

-6.4208 

-8.3393 

-8.4723 

5 . 3030 

3.8495 

-6.3225 

-6.4208 

3.8495 

2.9563 

Ordinal  Logistic  Regression:  Score  versus  logC.  log 

Link  Function:  Normit 

Response 

Informat  ion 

Variable 

Va  1  ue 

Count 

Score 

1 

1  1 

2 

9 

3 

18 

Total 

38 

Log i st ic 

Regression  Table 

Predictor 

Coef 

SE  Coef 

Z 

P 

Const ( 1 ) 

14.442 

4.039 

3.58 

0.000 

Const (2) 

15.339 

4.  106 

3.74 

0.000 

logC 

-9.825 

2.591 

-3.79 

0.000 

logT 

-7.028 

1 .894 

-3.71 

0.000 

Sex 

-0.3973 

0.2137 

1.86 

0.063 

LogT45sq 

1  .272 

1  .223 

1  .04 

0.299 

The  last  term  is  [Log(T/45)  ] A2 

Log- 1  ike  1 i hood  =  -30.139 

Test  that  all  slopes  are  zero:  G  =  19.821.  DF  =  4.  P  Value  =  0.001 


Goodness  of  - Fi t  Tests 


Method 
Pearson 
Dev i ance 


Chi -Square 
71.098 
60.278 


DF  P 

68  0.375 

68  0.736 


Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Pairs 

Number 

Percent 

Summary  Measures 

Concordant 

377 

82. 1% 

Somers ’  D 

0.64 

Discordant 

82 

1 7 . 9% 

Goodman- Kruskal  Gamma 

0.64 

Ties 

0 

0 . 0% 

Kenda 1 1 ' s  Tau -a 

0.42 

Total 

459 

100.0% 

MTB  >  print  m5 
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Data  Display 

Matrix  XPWX1 


16.3159 

16.5512 

10.2044 

-7.6127 

-0.1809 

0.6725 

16.5512 

16.8588 

10.3830 

-7.7411 

-0. 1885 

0.6955 

10.2044 

- 10.3830 

6.7116 

4.7776 

0. 1 195 

0.8497 

-7.6127 

-7.7411 

4.7776 

3.5871 

0.0832 

0.4464 

-0. 1809 

-0. 1885 

0. 1 195 

0.0832 

0.0457 

-0.0039 

0.6725 

0.6955 

-0.8497 

-0.4464 

-0.0039 

1 .4964 

MINITAB  PRINTOUTS  (WITHOUT  PIG  63) 


Model  LI 


Ordinal  Logistic  Regression:  Score  versus  logC,  Time,  sex 

Link  Function:  Norm  it 

Response  Information 

Variable  Value  Count 
Score  0  1 1 

1  9 

2  17 

Total  37 


Factor 

Informat 

ion 

Factor 

Leve 1 s 

Va 1 ues 

Time 

3 

10.  60. 

180 

sex 

2 

f  ema 1 e . 

male 

*  NOTE  *  37  cases  were  used 

*  NOTE  *  1  cases  contained  missing  values 

[DRS]  Female  Pig  63  is  an  outlier  and  was  dropped  from  this  analysis. 
Logistic  Regression  Table 


Predictor 

Coef 

SE  Coef 

Z 

P 

Const ( 1 ) 

13.8594 

3.66928 

3.78 

0.000 

Const (2) 

15.0153 

3.78178 

3.97 

0.000 

logC 

-  15.7332 

3.98037 

-3.95 

0.000 

Time 

60 

-8.64404 

2.34381 

3.69 

0.000 

180 

-  15. 1925 

3.89458 

-3.90 

0.000 

sex 

ma  1  e 

-  1 .81609 

0.922677 

-1.97 

0.049 

Time*sex 

60*male 

-0.632034 

1 . 13371 

-0.56 

0.577 

180*male 

1 .88576 

1 . 22680 

1 .54 

0.  124 

Tests  for 

terms  wi th 

more  than  1 

degree  of  freedom 

Term 

Ch i -Square 

DF  P 

Time 

15.2175 

2  0.000 
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Time*sex 


4.2498  2  0.119 


Log  L i ke 1 i hood  =  -24.252 

Test  that  all  slopes  are  zero:  G  =  30.070.  DF  =  6.  P-Value  =  0.000 

Goodness -ofFit  Tests 

Method  Chi  Square  DF  P 

Pearson  61.0914  64  0.580 

Deviance  48.5048  64  0.925 

Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


t 

Pa  i  rs 

Number 

Percent 

Summary  Measures 

Concordant 

390 

88.8 

Somers'  D 

0.78 

Discordant 

48 

10.9 

Goodman -Kruskal  Gamma 

0.78 

Ties 

Total 

1 

439 

0.2 

100.0 

Kendal  1 ' s  Tau-a 

0.51 

MTB  >  print  m2 

Data  Display 


Matrix  XPWX4 


13.4636 

13.8243 

-14.4182 

-8.22759 

-14. 1090 

-  1 .94651 

0.04925 

2.09104 

13.8243 

14.3018 

- 14.8744 

-8.48282 

-14.5525 

1 .99799 

0.03072 

2. 14581 

14.4182 

-  14.8744 

15.8434 

8.69882 

15. 1628 

1 .79405 

0.29641 

-1.95214 

8.2276 

-8.4828 

8 . 6988 

5.49343 

8.6392 

1 .29877 

-0.55390 

-  1 .38551 

14. 1090 

-14.5525 

15. 1628 

8.63922 

15. 1677 

2.03072 

-0.02941 

-2.52403 

1 .9465 

- 1 .9980 

1 .7940 

1 .29877 

2 . 0307 

0.85133 

-0.61463 

-0.86924 

0.0492 

0.0307 

0.2964 

-0.55390 

-0.0294 

-0.61463 

1 .28531 

0.61180 

2.0910 

2. 1458 

-1 .9521 

-  1 .38551 

-2.5240 

-0.86924 

0.61 180 

1 . 50503 

Model  L2 

Ordinal  Logistic  Regression:  Score  versus  logC,  Time,  sex 

Link  Function:  Norm  it 

Response  Information 

Variable  Value  Count 
Score  0  1 1 

1  9 

*  2  17 

Total  37 

+  Factor  Information 

Factor  Levels  Values 
Time  3  10.  60.  180 

sex  2  female,  male 

*  NOTE  *  37  cases  were  used 

*  NOTE  *  1  cases  contained  missing  values 

[DRS]  Female  Pig  63  is  an  outlier  and  was  dropped  from  this  analysis. 
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Logistic  Regression  Table 

Pred ictor 
Const ( 1 ) 

Const (2) 
logC 
Time 
60 
180 
sex 
male 

Tests  for 
Term  Chi-Square  DF  P 

Time  15.3604  2  0.000  * 

Log- Like  1 i hood  =  -26.484 

Test  that  all  slopes  are  zero:  G  =  25.607.  DF  =  4 .  P- Value  =  0.000 
Goodness - of -F i t  Tests 


Method 

Ch i -  Square 

DF 

P 

Pearson 

58.8779 

66 

0.721 

Dev i ance 

52.9675 

66 

0.877 

Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Coef 

SE  Coef 

Z 

P 

10.9275 

2.87594 

3.80 

0.000 

1 1 .9844 

2.97145 

4.03 

0.000 

- 12.6430 

3. 19098 

3.96 

0.000 

-7.34827 

1 .92722 

-3.81 

0.000 

-  1 1 .4296 

2.91644 

-3.92 

0.000 

-  1 . 19785 

0.488518 

2.45 

0.014 

terms  with 

more  than 

1  degi 

ree  of  freedom 

0.73 
0.73 
0.48 

MTB  >  print  m3 

Data  Display 

Matrix  XPWX5 


Pai  rs 

Concordant 

Discordant 

Ties 

Total 


Number 

379 

60 

0 

439 


Percent 

86.3 

13.7 

0.0 

100.0 


Summary 
Somers ’ 
Goodman 
Kendal  1 


Measures 

D 

Kruskal  Gamma 
s  Tau-a 


8.27103 

8.50165 

-9.0598 

-5.39319 

8.27102 

-0.675879 

8.50165 

8.82949 

-9.3671 

-5.57286 

-8.55020 

-0.705401 

9.05980 

-9.36708 

10. 1824 

5.91670 

9. 13973 

0.657970 

5.39319 

-5.57286 

5.9167 

3.71417 

5.45437 

0.354784 

8.27102 

-8.55020 

9.1397 

5.45437 

8.50561 

0.586322 

0.67588 

-0.70540 

0.6580 

0.35478 

0.58632 

0.238650 

Model  L3 


Ordinal  Logistic  Regression:  Score  versus  logC,  Time 

Link  Function:  Norm  it 
Response  Information 
Variable  Value  Count 


Score  0  1 1 

1  9 

2  17 

Total  37 
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Factor  Information 


Factor  Levels  Values 
Time  3  10,  60,  180 

*  NOTE  *  37  cases  were  used 

*  NOTE  *  1  cases  contained  missing  values 

[DRSJ  Female  Pig  63  is  an  outlier  and  was  dropped  from  this 


Logistic  Regression  Table 


Predictor 

Coef 

SE  Coef 

Z 

P 

Const ( 1 ) 

8.31350 

2.31076 

3.60 

0.000 

Const (2) 

9.21873 

2.38045 

3.87 

0.000 

logC 

-  10.2284 

2.65272 

-3.86 

0.000 

Time 

60 

-6.07721 

1 . 66820 

-3.64 

0.000 

180 

-9.31643 

2.44364 

-3.81 

0.000 

Tests  for 

terms  with 

more  than 

1  deg 

ree  of 

f  reedom 

Term  Chi 

-Square  DF 

P 

Time 

14.5493  2 

0.001 

Log-L i kel 

i hood  =  -  29 . 

.778 

Test  that 

all  si  opes 

are  zero: 

G  =  1 

9.019, 

DF  =  3. 

Goodness - of  - F i t  Tests 


Method 

Ch i -Square 

DF 

P 

Pearson 

77.7416 

67 

0.  174 

Deviance 

59.5555 

67 

0.729 

Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


Pai  rs 

Number 

Percent 

Summary 

Measures 

Concordant 

359 

81 .8 

Somers ' 

D 

0. 

64 

Discordant 

80 

18.2 

Goodman 

-Kruskal  Gamma 

0. 

64 

Ties 

0 

0.0 

Kendal  1 

’ s  Tau-a 

0. 

42 

Total 

439 

100.0 

MTB  >  print  m4 

Data  Display 


Matrix  XPWX6 


5.33961 

5.46658 

-6.04159 

-3.75441 

-5.56947 

5.46658 

5.66654 

-6.23213 

-3.86714 

-5.74398 

-6.04159 

-6.23213 

7.03695 

4.21734 

6.33423 

-3.75441 

-3.86714 

4.21734 

2.78289 

3.92894 

-5.56947 

5.74398 

6.33423 

3.92894 

5.97139 

Model  L4 

Ordinal  Logistic  Regression:  Score  versus  logC,  logT,  sex 

Link  Function:  Normit 


ana lys i s . 


0.000 


APPENDIX  B 


61 


Response  Information 


Variable  Value  Count 


Score  0  1 1 

1  9 

2  17 

Total  37 


Factor  Information 


Factor  Levels  Values 
sex  2  f ema 1 e .  male 

*  NOTE  *  37  cases  were  used 

*  NOTE  *  1  cases  contained  missing  values 


[DRS]  Female  Pig  63  is  an  outlier  and  was  dropped  from  this  analysis 


Logistic  Regression  Table 


Predictor 

Coef 

SE  Coef 

Z 

P 

Const ( 1 ) 

22.2598 

5.85398  3. 

80 

0.000 

Const (2) 

23.3509 

5.95337  3. 

92 

0.000 

logC 

13.5029 

3.42349  -3. 

94 

0.000 

logT 

- 10.2976 

2.68506  -3. 

84 

0.000 

sex 

male 

3.06349 

1.72433  -1. 

78 

0.076 

sex* logT 

ma  1  e 

1 .07283 

0.930951  1. 

15 

0.249 

Log-Likelihood  =  -25. 

921 

Test  that 

all  si  opes 

are  zero:  G  = 

=  26 

.733.  1 

Goodness  - 

of -Fit  Tests 

Method 

Ch i  Square 

DF  P 

Pearson 

55.4602 

66  0.819 

Dev i ance 

51.8415 

66  0.899 

Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


t 


% 


Pa  i  rs 

Number 

Percent 

Summary  Measures 

Concordant 

377 

85.9 

Somers'  D 

0.72 

Discordant 

62 

14.  1 

Goodman- Kruskal  Gamma 

0.72 

Ties 

0 

0.0 

Kendal  1  s 

Tau-a 

0.47 

Total 

439 

100.0 

MTB  >  print 

m5 

Data  Display 

Matrix  XPWX7 

34 . 2690 

34 . 8050 

-  19.7042 

-  15.6628 

-5.67531 

2.46103 

34.8050 

35.4426 

-20.0507 

-15.9283 

-5.75225 

2.48795 

- 19.7042 

20.0507 

1 1 . 7203 

8.9509 

2.70952 

■1 .09640 

-15.6628 

15.9283 

8 . 9509 

7 . 2096 

2.68421 

1 .21099 

-5.6753 

-5.7523 

2 . 7095 

2.6842 

2.97331 

1 .53557 

2.4610 

2.4879 

1 .0964 

-1.2110 

1 .53557 

0.86667 
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Model  L5 


Ordinal  Logistic  Regression:  Score  versus  logC,  logT,  sex 


Link  Function:  Norm  it 

Response  I  n forma t.  i  on 

Variable  Value  Count 
Score  0  1 1 

1  9 

2  17 
Total  37 

Factor  Information 

Factor  Levels  Values 
sex  2  female,  male 

*  NOTE  *  37  cases  were  used 


*  NOTE  * 

1  cases  contained  missing  va 

1  ues 

[DRS]  Female  Pig  63 

is  an  out  1 ier  and 

was  dropped 

from  this 

Log i st ic 

Regression 

Table 

Predictor 

Coef 

SE  Coef 

Z 

P 

Const ( 1 ) 

19.6521 

5.07896 

3.87 

0.000 

Const (2) 

20.6991 

5. 17367 

4.00 

0.000 

logC 

12.4093 

3. 12741 

-3.97 

0.000 

logT 

8.97914 

2 . 28809 

-3.92 

0.000 

sex 

ma  1  e 

1 .21073 

0.489222 

-2.47 

0.013 

Log  - L i ke 1 

i hood  =  -  26 

.640 

Test  that 

all  si  opes 

are  zero 

:  G  =  25 

.296.  DF  =  3. 

P-Value 

Goodness  - 

of-Fit  Tests 

Method 

Ch i -Square 

DF 

P 

Pearson 

60. 1095 

67  0.7 

12 

Dev i ance 

53.2794 

67  0.889 

Measures 

of  Association: 

(Between 

the  Response  Variable  and  Predicted  Probabilities) 

Pai  rs 

Number 

Percent 

Summary  Measures 

Concordant  374 

85.2 

Somers ' 

D 

0.70 

Discordant  65 

14.8 

Goodman - 

Kruskal  Gamma 

0.70 

Ties 

0 

0.0 

Kenda 1 1 ' 

s  Tau-a 

0.46 

Total 

439 

100.0 

MTB  >  print  ml 

Data  Display 


Matrix  XPWX3 


25.7958 

26.2335 

-15.6902 

-11.5781 

-  1 . 18443 

26.2335 

26.7668 

-15.9878 

-11.7947 

-  1.21315 

15. 6902 

15.9878 

9.7807 

7.0269 

0.68225 
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=  0.000 


-11.5781  -11.7947  7.0269  5.2353  0.48445 

-1.1844  -1.2131  0.6823  0.4845  0.23934 

Model  L6 


Ordinal  Logistic  Regression:  Score  versus  logC,  logT 

Link  Function:  Normit 

Response  Information 

Variable  Value  Count 
Score  0  1 1 

1  9 

2  17 

Total  37 

*  NOTE  *  37  cases  were  used 

*  NOTE  *  1  cases  contained  missing  values 

[DRS]  Female  Pig  63  is  an  outlier  and  was  dropped  from  this  analysis. 

Log i s t i c  Regress  ion  Table 

Predictor 
Const ( 1 ) 

Const (2) 
logC 
logT 

^  i  ke  1  i  hood  =  -30.001 

2.  P- Value  =  0.000 


Coef 

SE  Coef 

Z 

P 

15. 1659 

4.05345 

3.74 

0.000 

16.0604 

4. 12562 

3.89 

0.000 

-9.88150 

2.54186 

3.89 

0.000 

-7.22086 

1 . 88084 

3.84 

0.000 

i  hood  =  -  30 . 
al  1  slopes 

.001 

are  zero: 

G  =  1 

18.572. 

Goodness- of -Fi t  Tests 

Method  Chi-Square  DF  P 
Pearson  79.7589  68  0.156 
Deviance  60.0028  68  0.744 


Measures  of  Association: 

(Between  the  Response  Variable  and  Predicted  Probabilities) 


r 

* 


Pa  i  rs 

Number 

Percent 

Summary  Measures 

Concordant 

352 

80.2 

Somers'  D 

0.61 

Discordant 

86 

19.6 

Goodman -Kruskal  Gamma 

0.61 

Ties 

1 

0.2 

Kendal  1 ' s  Tau-a 

0.40 

Total 

439 

100.0 

MTB  >  print 

m6 

Data  Display 

Matrix  XPWX8 

16.4305 

16. 

6894 

-  10. 1366 

-7.59962 

16.6894 

17. 

0207 

-  10.3226 

7.73621 

-10. 1366 

10. 

3226 

6.4611 

4.66624 

-7.5996 

-7. 

7362 

4 . 6662 

3.53758 
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